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This  vork  ves  performed  by  the  UAC  Research  laboratories  for 
the  Aeronautical  Syst—  Division,  Vtight-Patterson  Air  Foret  But, 
under  Contract  A7  3$(6l6)-6800,  Project  Bo.  1366,  Task  Bo.  13660$. 

Tfc it  Project  and  Task  It  a  part  of  the  Air  Foret  Syttaat  Co— nil  'a 
Applied  Autarch  Program  75OA,  Itechanics  of  Flight/^  This  contract 
vai  iasued  to  the  UAC  Research  Laboratories  on  the  "basis  of  their 
response  to  Purchase  Bequest  Ho.  65496  and  vork  vu  initiated  under 
the  contract  on  April  4,  i960 . 


Included  among  those  who  co-operated  in  the  research  and 
preparation  of  the  report  vere:  Mr.  M.  Schweiger,  Head,  Propulsion 
Section,  Mr.  C.  E.  Kepler,  Supervisor,  Inlet  Group,  and  Mr.  R.  L. 
O'Brien,  Project  Engineer.  The  vork  vas  monitored  by  Mr.  P.  H. 
Kutschenreuter  of  the  Flight  Dynamics  laboratory,  ASD. 


The  characteristics  of  supersonic  tuibulent  boundary  layers  on  cooled 
surfaces  haring  adverse  pressure  gradients  were  investigated  to  provide 
design  information  for  hypersonic  inlets.  Tests  were  conducted  at  Kach 
mafoers  of  3  and  6  with  values  of  vail  temperature  as  low  as  300  of  the 
adiabatic  recovery  temperature .  The  adverse  pressure  gradients  were  gen¬ 
erated  by  two-dimensional  models  including  a  circular-arc  surface,  an 
isentropic-compression  surface,  and  an  oblique  shock  generator.  Velocity 
and  temperature  profiles  through  the  boundary  layer  were  obtained  and  the 
integral  boundary  layer  parameters  were  evaluated  along  the  surface.  An 
analytical  investigation  was  also  conducted  to  provide  a  basis  for  evalu¬ 
ating  the  experimental  results. 

It  was  found  that  cooled,  turbulent  boundary  layers  in  adverse  pres¬ 
sure  gradients  are  thinner  and  less  distorted  than  uncooled  boundary  layers, 
that  the  ratio  of  Stanton  nuaber  to  skin-friction  coefficient  is  higher  than 
the  flat-plate  value,  and  that  the  total  amount  of  heat  removed  is  greater 
than  that  predicted  by  integrating  local  heat-transfer  rates  based  on  flat- 
plate  conditions  along  the  surface. 
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X.  Amps  of  ft  supersogie  tuibulsat  bousdar?  3*tr  in  an  adverse 
(positive)  pressure  gradient  ornate*  by  a  curved  compression  surface  li 
laproved  by  the  application  of  mil  coolli*  tad  a  significant  reduction  la 
boundary  layer  thickness  la  bbtainsd.  Calculation  of  hypersonic  inlet 
parfozmaoa  by  aathoda  vhich  mold  proparly  account  for  thaaa  changes  la 
ilia  boundary  layar  characteristics  vith  vail  cooling  (for  exsmple,  tha 
aethod  of  Ref.  2)  would  Indleata  that  vail  coollag  will  improve  tha  total 
pressure  recovery  of  hyparsoale  Inlet  a.  The  lncreaae  la  praasura  re  core  ry 
la  a  result  of  two  factors j  tha  boundary  layar  vith  vail  cooling  can  nego¬ 
tiate  a  higher  pressure  rise  without  separating,  and  the  reduction  of 
boundary  layer  thickness  vith  cooling  allows  higher  maximum  inlet  contrac¬ 
tion  ratios  to  be  Obtained. 


2.  The  ratio  of  Stanton  nuriber  to  skin-friction  coefficient  for  a 
supersonic  turbulent  boundary  layer  on  a  cooled  flat  plate  is  independent 
of  Mach  number  in  the  range  froai  3  to  6.  For  the  configuration  tested  (a 
flat  plate  downstrean  of  a  streanvise  step  in  vail  temperature  distribution), 
this  ratio  is  higher  than  that  predicted  on  the  basis  of  the  Reynolds -Colburn 
analogy. 


3.  The  ratio  of  Stanton  number  to  skin-friction  coefficient  for  a 
supersonic  .turbulent  boundary  layer  in  an  adverse  pressure  gradient  is  lar¬ 
ger  that  on  a  flat  plate  and  the  deviation  frost  the  flat-plate  value 
Increases  with  the  magnitude  of  the  gradient. 

4.  The  total  amount  of  heat  removed  frost  a  supersonic  turbulent 
boundary  layer  on  a  cooled  surface  having  an  adverse  pressure  gradient  is 
greater  than  that  predicted  by  integrating  local  heat-transfer  rates  based 
on  flat-plate  conditions  along  the  surface. 


Manuscript  released  by  the  authors  November  1961  for  publication  as  an  A8D 
Technical  Report 


1.  Adlltlomal  investigations  of  supsrsonie  tattoo loot  boundary  Upon 
on  ooolad  smrlhots  having  advtrs*  pressure  gradients  should  bs  conducted 
to  provide  sort  detail**  experimental  dot*  for  om  in  the  development  of 
generalised  procedures  for  predicting  boundary  lsyer  characteristics  la 
hypersonic  Inlet* ,  since  the  results  of  this  study  (see  CWCUJ8I0*  l)  in¬ 
dicate  th at  vail  cooling  should  bore  a  significant  favorable  effect  on 
the  performance  of  hypersonic  inlets.  Specifically,  the  following  charac¬ 
teristics  should  be  Investigated: 

a.  The  effect  of  Reynolds  matter  on  the  shape  parameter 
for  boundary  layers  on  cooled  vails. 

b.  Ifce  effect  of  wall  cooling  on  boundary  layer  separation 
criteria. 


c.  The  effect  of  vail  cooling  on  duct  flows  comprised 
entirely  of  boundary  layer. 


2.  Investigations  Should  be  conducted  to  determine  the  characteristics 
of  supersonic  laminar  boundary  layers  on  cooled  vails  having  pressure  grad¬ 
ients  which  simulate  those  of  hypersonic  inlets  designed  for  supersonic- 
burning  ramjets. 


3.  Improved  experimental  techniques  should  be  developed  for  measuring 
boundary  layer  total  temperature  profiles  adjacent  to  cooled  vails  so  that 
more  accurate  data  can  be  obtained  in  future  studies.  Techniques  for  meas¬ 
uring  vail  temperature  differences  accurately  should  also  be  developed  so 
that  direct  measurements  of  wall  heat-transfer  rates  can  be  obtained. 


4.  Future  hypersonic  inlet  investigations  should  be  conducted  with 
vail  cooling  to  provide  proper  simulation  of  the  boundary  layer  character¬ 
istics  in  such  inlets. 


ft*  toaster  l«r»r  aevalopsnt  os  the  ualle  of  supersonic  and  typer- 
•oalc  Inlet*  feu  a  large  influence  oa  the  performance  of  Maoh  amber, 
air-breathing  propulsion  systems.  Ill*  boundary  layer  characteristics  Halt 
th*  amount  of  supersonic  ooupressloe  and,  for  systems  suplcyii*  subsonic 
eoAastloa,  lnflaeno#  th*  pressure  recovery  through  th*  norml  shock.  Al¬ 
though  detailed  boundary  layer  growth  characteristics  haviig  direct  appli¬ 
cation  to  the  design  of  supersonic  inlets  with  uncooled  wall*  were  Obtained 
in  a  previous  investigation  at  th*  QIC  Research  Laboratories  (Ref.  1),  these 
data  are  not  generally  applloSbl*  to  hypersonic  inlets  where  vail  coolly 
vlll  be  required.  Recent  studies  of  hypersonic  inlets  jtt  me  (Ref.  2)  have 
shown  that  wall  cooling  has  a  significant  effect  on  th*  performance  of  hyper 
sonic  inlets.  This  is  to  be  expected  since  boundary  layer  growth  and  sep¬ 
aration  fhenouena  rlajr  a  large  role  in  determining  th*  perforaance  of 
hypersonic  inlets.  Xn  such  inlets,  the  boundary  layer  constitutes  the 
■ajar  portion  of  the  flow  at  the  throat. 


The  objective  of  the  investigation  described  herein  was  to  Obtain 
detailed  information  on  supersonic  turbulent  boundary  layer  growth  over 
cooled  walla  in  adverse  pressure  gradients  to  provide  a  note  rational  basis 
for  the  design  of  hypersonic  inlets. 
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The  only  available  experimental  data  on  nqperaondo  turbulent  boundary 
layer  characteristics  over  cooled  or  heated  valla  In  adrerae  pressure  grad¬ 
ients  la  contained  In  Ref.  3*  This  work  was  Halted  to  a  single  case  of  a 
strong  oblique  shock  Interaction  on  a  vail  which  was  only  moderately  cooled 
or  heated  (0.93<TW/T-  <1.26),  and  the  data  consisted  primarily  of  vail 
static  pressure  distributions.  From  this  data  the  author  deduced  that  a 
small  Increase  of  the  separation  pressure  ratio  vas  obtained  with  cooling. 

A  prior  survey  of  the  literature  (Ref.  4)  contains  data  for  subsonic  cooled 
turbulent  boundary  layers  In  adverse  pressure  gradients  but  no  comparable 
results  for  supersonic  flov.  Because  of  this  lack  of  experimental  data, 
theoretical  calculation  procedures  for  cooled  supersonic  boundary  layers 
In  adverse  pressure  gradients,  such  as  that  of  Ref.  5,  are  extensions  of 
the  procedures  developed  for  Incompressible  flov,  e.g.,  Ref.  6. 

Limited  data  exist  for  turbulent  boundary  layers  on  cooled  flat  plates 
(zero  pressure  gradient)  at  high  supersonic  Mach  nuabers,  see  Refs.  7 
through  10.  These  data  are  not  in  complete  agreement  vlth  existing  theories 
(see,  for  example,  Refs.  5>  11,  and  12);  the  data  indicate  a  trend  of  de¬ 
creasing  skin  friction  with  vail  cooling,  vhereas  the  theories  all  predict 
an  Increase. 

Supersonic  turbulent  boundary  layer  growth  over  insulated  walls  In 
adverse  pressure  gradients  has  been  investigated  quite  extensively  for  the 
case  of  the  severe  adverse  pressure  gradients  caused  by  oblique  or  normal 
shocks  (see,  for  example,  Refs.  13  through  23).  For  this  class  of  gradients, 
the  mixing  of  flow  from  the  free  stream  into  the  boundary  layer  has  little 
effect  on  the  shape  of  the  boundary  layer  over  the  short  length  of  pressure 
rise,  and  several  theories  (Refs.  2k  through  28)  have  been  developed  to 
determine  the  maximum  pressure  rise  that  can  be  obtained  without  separation. 
The  only  known  data  on  supersonic  boundary  layers  in  adverse  pressure  grad¬ 
ients  not  created  by  shocks  are  contained  in  Refs.  1,  29,  and  30.  The  theory 
presented  in  Ref.  5  is  applicable  for  very  mild  pressure  gradients  of  this 
type;  however,  for  large  gradients,  the  theory  severely  underestimates  the 
pressure  rise  which  can  be  obtained  without  separation  of  the  boundary  layer. 
This  discrepancy  between  theory  and  experiment  (see  Ref.  1  for  detailed  dis¬ 
cussion)  occurs  because  the  theory  does  not  account  for  increased  mixing  of 
the  freestream  flow  into  the  boundary  layer.  The  effect  of  this  increased 
mixing  is  included  in  the  "lag-length  procedure"  developed  in  Ref.  1  for 
calculating  supersonic  turbulent  boundary  layer  growth  in  adverse  pressure 
gradients.  This  theory,  which  is  based  on  experimental  data  obtained  from 
tests  in  which  the  wall  temperature  was  in  equilibrium  for  zero  heat  trans¬ 
fer  (adiabatic  wall),  is  modified  in  the  following  section  to  account  for 
wall  cooling. 
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AHALmCAL  aVESTIOATIOlf 


Prediction  of  Boundary  Layer  Growth  and 
Beat  Tranafer  on  a  Plat  Plate 


The  theoretical  aethoda  for  predicting  turbulent  boundary  layer  growth 
and  heat-transfer  charaeteristiea  on  a  flat  plate  employed  in  this  report 
were  derived  from  the  akin-friction  coefficient  relationship  obtained  from 
Ref.  5.  In  this  reference  the  incompressible  Ludwieg-ffillmann  skin-friction 
equation  (Ref.  31)  is  modified  for  compressibility  by  using  the  Eckert  ref¬ 
erence  enthalpy  method  of  Ref.  32.  The  relationship  for  skin-friction  co¬ 
efficient  used  in  Ref.  5  can  be  rewritten  in  the  following  forms 


C 


f 


0.246eH56lHi 


(1) 


where  Hj  is  the  inconpressible  shape  parameter  and  Tref  is  the  reference 
temperature  as  defined  in  Ref.  32.  In  this  derivation  the  power  law  for 
viscosity,  ^~T  •  ,  has  been  used.  After  integrating  the  von  Karman 

boundary  layer  momentum  equation  (employing  Eq.  1  to  evaluate  Cf  ),  the 
boundary  layer  momentum  thickness, 9  ,  may  be  evaluated  by  the  following 
expression: 


(Po»o  \'0,212 

/  Jso  f  628 

\  Po  1 

'tref' 

(x)' 


0.788 


(2) 


where  the  boundary  conditions  are  9-0  at  x=0  . 
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T!»  boundary  layer  displacement  thickness,^  ,  vu  evaluated  from  the 
theoretical  value  of  8  and  the  theoretical  value  of  H  *8 */0  for  a  cooled 
1/7 -power  profile  ae  tabulated  in  Ref.  33*  This  reference  employed  a  mod¬ 
ification  of  the  Crocco  energy  Integral  to  relate  total  temperature  and 
velocity  for  evaluating  the  Integral  boundary  layer  parameters . 

The  equation  relating  heat -transfer  rate,  q  ,  to  the  total  beat  removed, 
0,  Is 


q 


dQ 

dx 


s,/>ouocp  {tr-V 


(3) 


The  total  heat  removed,  Q  ,  is  related  to  the  total  enthalpy  thickness,  6 
by  the  following  relation 


0  ■  />ou ohT0^ 


where  <f>  is  defined  by 


(5) 


Thus,  for  a  flat  plate,  <f>  is  directly  proportional  to  the  total  heat  re¬ 
moved  from  the  boundary  layer.  If  the  ratio  of  the  Stanton  number  (heat- 
transfer  coefficient)  to  skin-friction  coefficient,  St/Cf  ,  is  a  constant, 
Eq.  3  may  be  integrated  directly  to  evaluate  the  total  enthalpy  thickness 
as  a  function  of  distance  along  the  surface.  The  resulting  expression  is 

■*”*°+2lr  Ur  -jr-h-Bo)  <6> 

1  70  To 
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vbtrt  tb>  boundary  conditions  ar*  »t  8  *  90  .  The  ratio  of  Stanton 
aurtwr  to  ikln-frletloc  coefficient  given  by  the  Reynolds -Colburn  analogy 
is  equal  to  0(635  for  a  Prendtl  nuaber  of  0.70. 
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Prediction  of  Boundary  Layer  Growth  and  Seat 
in  an  Adverse  Pressure  Gradient 


Transfer 


The  asthod  used  in  this  report  to  predict  boundary  l»v*r  4 

adverse  pressure  gradient  is  based  on  the  von  “ 

equation  for  compressible  flow  (Ref.  34) »  aomentum  integral 


d  Q  * 


+ 


H 


"w  / 


(7) 


TZ? 1*2  0?  tM*.”11  »Tlc  »—»•  •* 

^re  It  i.  concluded  that  the  nuLnL  tSkn^e  boL^ 

«  frsr ls1 

vcrs^r5”'*’1”  ■urfaee  1,1  sisns  2, 

suined  to  Increase  at  the  flat -plate  value  over  the  1s  as“ 

pre88r 8radlent'  £“££  *£. 

rate  of  mass  flow  into  the  boundary  layer  increases  rsr^w  gD1®ngth>  the 

rr0 18  a 

«-  SJSSfiST rT  V: llvlded  * 

G=  H/H',/7  ras  defied;  the  ™i£Af  *,/.  parMKter' 

For  an  uncooled  boundary  laver  a  a  nr  n^+J7/  -,yere  obtained  from  Ref.  35. 

a  value  of  G>l  implies  a  distorted  profile.  *  proflle  and 

surfaced  !5£'f  iTT  ^  “  «W 

for  the  Insula ted-valAaee  uslnTthe  la.  lenrtA  ^  6  ras  flrst  oraiPuted 
Ref.  1.  The  effect  of  ooo!l£  „  S  ™  t£f LSfJV8  4erelop'd  ltt 

ratl°  °f  -  SXTSSX  SS5S  S 
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gcooled 


^ADIABATIC 


Gc,WVEP  qAT. . £h* IS,. 

^ADIABATIC  FLAT  PLATE 


) 


(8) 


Eero  the  a ■  sumption  la  made  that  the  ratio  of  the  shape  parameter  on  the 
cooled  surface  to  the  shape  parameter  on  the  adiabatic  (insulated)  surface 
Is  the  same  for  the  distorted  profile  (resulting  from  the  adverse  pressure 
gradient)  as  it  la  for  a  flat  plate  profile.  For  the  theoretical  shape 
parameters  tabulated  In  Ref.  33  for  various  power  lav  profiles,  calculations 
of  the  ratio  of  the  shape  parameter  on  a  cooled  surface  to  the  shape  para¬ 
meter  on  an  adiabatic  surface  at  given  values  of  Mach  number  and  wall  temper¬ 
ature  indicate  that  the  value  of  this  ratio  Is  approximately  constant  for 
power-lav  exponents  ranging  from  l/5  to  l/ll.  This  assumption  is  Investi¬ 
gated  further  in  APPENDIX  III  where  an  attempt  Is  made  to  evaluate  the 
validity  of  Eq.  8  by  an  analysis  of  the  data  obtained  In  this  Investigation. 

The  value  of  G  for  a  cooled  flat  plate  vas  obtained  from  the  experimental 
data  for  the  flat-plate  studies  described  in  a  later  section  of  this  report. 

The  skin-friction  coefficient  was  computed  from  Eq.  1  which  Includes 
the  effects  of  a  alstorted  profile  In  the  incompressible  shape  parameter, 

H|  .  The  momentum  thickness  vas  then  computed  by  a  step-by-step  Integra¬ 
tion  of  Eq.  7  and  the  displacement  thickness  vas  computed  from  Q  using  the 
value  of  G  from  Eq.  8. 

Two  reference  values  of  heat  transfer  were  computed  for  comparison  with 
the  experimental  data.  The  first  value,  herein  termed  the  "lover  probable 
limit”  of  total  heat  transfer,  was  computed  by  assuming  that  the  local  heat- 
transfer  rate  at  each  point  on  the  surface  is  the  same  as  that  for  a  flat  - 
plate  profile  at  the  same  wall  temperature,  Mach  number,  and  Reynolds  number 
based  on  momentum  thickness.  This  is  the  method  suggested  In  Ref.  5  and  is 
commonly  used  in  boundary  layer  and  heat  transfer  calculation  procedures.  The 
second  reference  value,  herein  termed  the  "upper  probable  limit"  of  total 
heat  transfer,  was  cooputed  by  assuming  that  the  total  heat  removed  up  to  a 
given  station  In  an  adverse  pressure  gradient  is  equal  to  the  value  that 
would  have  been  removed  from  the  flat -plate  profile  having  the  same  9  , 
assuming  that  the  boundary  layer  developed  over  a  flat  plate  having  constant 
values  of  the  wall  temperature  and  Mach  number  equal  to  those  at  the  given 
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station.  In  all  eases  calculations  show  that  this  second  value  of 
total  heat  removed  is  much  larger  than  the  first  value.  Thus  the 
second  method  (higher  total  beat  removed)  would  result  in  a  lover 
average  total  temperature  of  the  air  within  the  boundary  layer  than 
the  first  method  (lover  total  heat  removed).  If  it  16  assumed  that 
the  rate  at  which  heat  is  transferred  from  a  boundary  layer  to  a 
wall  is  a  function  of  the  difference  between  the  average  temperature 
of  the  air  within  the  boundary  layer  and  the  wall  temperature,  then 
the  heat-transfer  rate  for  the  actual  boundary  layer  must  be  greater 
than  the  value  used  in  evaluating  the  lover  probable  limit  (first 
method)  and  less  than  the  value  required  to  establish  the  upper 
probable  limit  (second  method). 


The  total  heat  removed  at  any  station  along  the  surface  was  ,  , 
computed  in  terms  of  the  corrected  total  enthalpy  thickness , <£c  *  *  <p 

using  the  following  equations:  0  0 


For  the  lower  probable  limit, 


/  Jr. 
Wo  \Tt0 


(9) 


where  St  is  evaluated  for  the  flat-plate  value  as  determined  by  Reg  , 
Mw  ,  and  Vtt0  • 

For  the  upper  probable  limit,  s  4>'c^ 


^Uw 

?ouo 


'-So) 


(10) 
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vhere  St/Cf  *  0.5  Pr“2^5  s  0.635. 

The  method  for  predicting  the  turbulent  boundary  layer  growth  and 
heat -transfer  characteristics  developed  in  this  section  were  employed 
to  calculate  the  theoretical  values  presented  with  the  data  in  subse¬ 
quent  sections. 

The  method  for  calculating  turbulent  boundary  layer  growth  over 
a  cooled  compression  surface  presented  in  Ref.  5  was  also  employed  for 
the  case  of  the  flow  over  an  isentropic  compression  surface  at  Mach  6 
(one  of  the  configurations  tested).  The  results  indicate  that  the 
boundary  layer  would  separate  at  a  total  turning  angle  of  18  deg, 
whereas  the  method  described  above  did  not  predict  separation  at  any 
point  on  the  surface  having  a  total  turning  angle  of  34  deg.  Schlieren 
photographs  of  the  flow-over  model,  presented  in  Fig.  26,  give  no 
indication  of  boundary  layer  separation  at  any  point  on  the  surface  of 
the  model.  A  brief  description  of  the  lag-length  method  which  was 
employed  for  predicting  boundary  layer  separation  and  a  comparison  of 
this  method  with  Ref.  5,  has  been  presented  in  the  RESUME  OF  PREVIOUS 
WORK. 


EXPERIMENTAL  INVESTIGATION 


Description  of  Test  Equipment 

Tunnel 

Tests  to  determine  turbulent  boundary  layer  growth  over  cooled 
walls  in  adverse  pressure  gradients  were  conducted  in  the  UAC  Research 
Laboratories  4.5-x  3.5-in.  continuous-flow  supersonic  wind  tunnel.  Dry 
air  is  supplied  to  this  tunnel  at  a  total  temperature  of  approximately 
400  F  and  at  stagnation  pressures  up  to  400  psig  and  is  exhausted  to  , 
vacuum  pumps.  The  tunnel  is  equipped  with  interchangeable  single  super¬ 
sonic  nozzle  blocks  (Mach  3  or  Mach  6)  and  a  centerline  block  on  which 
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the  boundary  layer  approaching  the  aodela  develops  naturally.  The 
Reynolds  nuabers  in  the  tunnel  test  section  were  O.32  x  10^  per  in. 
at  Mach  3  and  O.59  x  106  per  in.  at  Mach  6. 


Models 


Three  models  were  employed  for  the  tests:  a  flat  plate,  a  compres¬ 
sion  surface  having  a  constant  radius  of  curvature  and  a  total  turning 

angle  of  26  deg  (the  circular-arc  model),  and  a  Mach  6  isentropic  compres¬ 
sion  surface  having  a  total  turning  angle  of  3k  deg  (the  isentropic- 
surface  model).  The  inviscid  isentropic  compression  contour  was  corrected 
for  boundary  layer  displacement  effects  using  the  method  outlined  in  the 
ANALYTICAL  INVESTIGATION  to  arrive  at  the  theoretical  displacement  thick¬ 
ness.  For  one  series  of  tests,  a  shock  generator  which  could  be  oriented 
at  various  angles  of  attack  (0<  a <10  deg)  with  respect  to  the  free  stream 
was  employed  to  create  a  shock-reflection  interaction  on  the  flat  plate. 
Photographs  of  the  models  installed  in  the  tunnel  are  presented  in  Figs  1 
2,  and  3.  The  model  surfaces  were  fabricated  from  copper  to  provide  good  ' 
seat-transfer  characteristics  between  the  measured  boundary  layer  on  the 
upper  surface  of  the  models  and  the  coolant  (liquid  »2)  on  the  lower  sur- 

T?e  "?Cl?.ye”  a‘t“ctled  t0  “  triple-pass  heat  exchanger  as  shown 
m  Fig.  4.  The  liquid  nitrogen  from  a  throttled  125  psig  supply  was  first 
passed  between  the  copper  models  and  the  closely  spaced  inner  plates  to 
create  high  coolant  velocities;  then  the  coolant  flow  made  a  double  pass 
beneath  the  models  to  minimize  heat  conduction  from  the  hot  walls  of  the 
tunnel  to  the  cold  model.  The  spacing  between  the  copper  models  and  the 
inner  plate  and  the  thickness  distributions  of  the  copper  models  were 
designed  so  that  a  constant  wall  temperature  could  be  maintained  on  the 
air  side  of  the  models.  The  heat- transfer  rates  employed  in  these  cal¬ 
culations  were  based  on  the  lower  probable  limit  of  heat  transfer  as 
defined  in  the  ANALYTICAL  INVESTIGATION. 


The  sides  of  the  compression-surface  models  were  cut  back  from  the 
walls  of  the  tunnel  in  order  to  minimize  the  effect  of  tunnel  side-wall 
boundary  layer  on  the  flow  over  the  models.  This  design  was  found  to 
be  necessary  for  similar  tests  reported  in  Ref.  1.  The  models  were 
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Instrumented  at  various  stations  with  both  surface  thermocouples  and 
static  pressure  orifices  for  measuring  vail  temperature  and  static 
pressure  distributions. 

Probe 

The  probe  employed  to  measure  the  boundary  layer  (see  Pigs.  1, 

2,  or  3)  vas  driven  by  a  mechanism  vhich  allowed  the  probe  to  traverse 
at  any  angle  In  the  longitudinal  vertical  plane  so  that  the  survey 
could  be  made  normal  to  the  model  surface.  The  probe  position  vas 
measured  to  an  accuracy  of  0.002  in.  The  probe  head  permitted  simul¬ 
taneous  measurement  of  pitot  pressure  and  temperature.  A  sketch  of 
the  probe  Is  presented  in  Pig.  3.  An  iron-constantan  thermocouple 
junction  vas  located  Inside  the  probe  head  (vhich  had  an  opening 
0.010  in.  high),  and  a  small  fraction  of  the  flow  vas  bled  through 
the  probe  to  provide  airflow  past  the  thermocouple  Junction.  The 
method  of  evaluating  the  local  total  temperature  from  the  measured 
temperature  is  discussed 'in  APPENDIX  I. 

The  effect  of  the  small  amount  of  bleed  flow  through  the  probe  on 
the  measured  pitot  pressure  vas  evaluated  on  the  basis  of  normal  shock 
inlet  data.  Typical  pressure-recovery,  weight- flow  data  for  a  normal 
shock  inlet  at  Mach  numbers  of  1.6  and  2.5  are  shown  in  Pig.  5.  The 
size  of  the  bleed  hole  in  the  probe  vas  selected  to  pass  a  0f 

15  to  200  of  the  ideal  weight  flow  through  the  probe.  Over  this  range 
of  weight  flows,  the  measured  pressure  recovery  of  a  normal  shock  inlet 
is  identical  to  the  dead-ended  pitot  pressure.  An  experimental  verifi¬ 
cation  of  this  characteristic  was  performed  by  surveying  the  boundary 
layer  on  a  flat  plate  with  both  a  conventional  pitot  probe  and  the 
combination  pressure- temperature  probe.  The  measured  pressures  were 
identical. 


Data  Reduction  Procedure 

Conventional  data  reduction  procedures  employed  in  the  evaluation 
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of  boundary  layer  parameters  on  flat  platae  cannot  be  used  for 
boundary  layers  on  curved  surfaces  because  of  the  large  static 
pressure  gradients  norml  to  the  surface.  These  gradients  provide 
the  centripetal  force  necessary  to  turn  the  flow.  Although  the 
wall  static  pressure  was  measured,  the  static  pressure  distribution 
through  the  boundary  layer  could  not  be  measured  accurately  by  a 
static  pressure  probe  because  of  flow  angularity.  Consequently,  the 
method  of  estimating  the  static  pressure  distribution  developed  in 
Kef.  1  was  employed.  The  essential  features  of  this  method  are  as 
follows:  A  typical  pitot  pressure  profile  normal  to  a  curved  surface 
increases  as  the  distance  from  the  wall,  y  ,  increases,  reaches  a 
maximum,  and  then  decreases  as  the  probe  traverses  into  the  free  stream 
where  there  is  also  a  static  pressure  and  Mach  number  gradient.  The 
edge  of  the  boundary  layer  is  usually  10  to  200  beyond  the  y  value 
for  the  peak  pitot  pressure.  The  static  pressures  Just  outside  the 
boundary  layer  were  evaluated  from  the  measured  pitot  and  free-stream 
total  pressures.  The  shape  of  the  static  pressure  distribution  through 
the  boundary  layer  was  obtained  by  fairing  a  curve  which  matched  the 
value  and  slope  of  the  free-stream  static  pressure  distribution  and  the 
value  and  slope  at  the  wall.  The  slope  at  the  wall  was  assumed  zero 
on  the  basis  of  a  zero  velocity  and  the  point  in  the  free  stream  at 
which  the  static  pressure  curves  were  matched  was  usually  20-300  beyond 
the  y  value  for  the  peak  pitot  pressure.  The  estimated  static 
pressure  variation,  the  measured  pitot  pressures,  and  the  measured  total 
temperature  could  then  be  employed  to  determine  all  the  other  local  flow 
properties . 

Three  integral  boundary  layer  parameters  were  evaluated  to  describe 
the  boundary  layer  characteristics:  the  mass  parameter  or  displacement 
thickness,  8*  ;  the  momentum  parameter  or  momentum  thickness,  9  ;  and 
the  energy  parameter  or  total  enthalpy  thickness,  4>  .  The  first  two  are 
the  conventional  parameters;  the  third  is  a  measure  of  the  total  heat 
removed  by  cooling  prior  to  the  station  at  which  the  profile  is  evaluated. 
This  parameter  was  employed  to  estimate  the  heat-transfer  rates  since 
measurements  of  the  temperature  gradient  across  the  cooled  copper  wall  and 
heat  rise  in  the  coolant  could  not  be  determined  accurately.  The  method 
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of  evaluating  these  integrals  is  presented  pictorially  in  Fig.  6. 

In  all  cases  the  parameters  are  based  on  the  ideal  conditions 
(primed  quantities)  that  would  exist  at  the  vail  (subscripted  w  ) 
as  determined  from  the  measured  vail  static  pressure,  free-stream 
total  pressure,  and  free-stream  total  temperature .  These  ideal  vail 
conditions  were  used  rather  than  the  conditions  at  the  edge  of  the 
boundary  layer  because  the  large  variation  of  the  conditions  outside 
the  boundary  layer  do  not  provide  a  convenient  set  of  reference  condi¬ 
tions.  It  can  be  seen  from  Fig.  6  that  8”  and  6  as  defined  in  this 
report  reduce  to  the  standard  definitions  of  these  quantities  for  the 
case  of  a  zero  pressure  gradient  normal  to  the  surface.  For  the  cases 
considered  herein  the  existence  of  a  normal  pressure  gradient  introduces 
additional  terms  into  the  von  Kdrmin  momentum  integral  equation.  Due  to 
the  complexity  of  these  terms  and  the  uncertainty  introduced  into  the 
data  reduction  procedure  by  the  necessity  of  estimating  the  static 
pressure  distribution  through  the  boundary  layer,  the  magnitude  of  this 
effect  is  difficult  to  evaluate.  ThiB  term  was  neglected  on  the  basis 
of  the  conclusions  of  Ref.  1  where  it  was  shown  for  a  sizeable  data 
sample  that  the  von  K&rmAn  momentum  integral  equation  in  the  form  of 
Eq.  7  predicted  the  boundary  layer  growth  characteristics  on  curved 
surfaces.  The  magnitude  of  these  terms  is  zero  for  the  case  of  a  zero 
normal  pressure  gradient.  The  definition  of  <p  ,  however,  allows  an 
exact  calculation  of  the  total  heat  removed  (using  Eq.  9),  when  St 
is  known  exactly. 


Discussion  of  Experimental  Results 


Mach  3  Studies 
Flat_Plate_ 

Tests  were  conducted  at  Mach  3  with  the  flat-plate  model  installed 
to  establish  the  initial  flow  conditions  for  subsequent  tests.  Three 
coolant  rates  were  employed:  zero,  the  maximum,  and  one  intermediate 
value.  The  conditions  which  correspond  to  these  three  coolant  flow 
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rates  are  hereafter  referred  to  as  "uncooled",  "maximum  cooling", 
and  "intermediate  cooling".  The  max  1  sun  coolant  flow  rate  was 
established  as  the  rate  for  which  liquid  nitrogen  was  exhausted  from 
the  heat  exchanger;  increasing  the  coolant  flow  rate  beyond  this 
value  provided  little  decrease  in  the  wall  temperature  ratio.  Figure  7 
shows  the  distributions  of  the  wall  Mach  number  and  wall  temperature 
for  each  of  the  three  coolant  conditions.  It  can  be  seen  that  the 
Mach  numbers  and  wall  temperatures  were  uniform  over  the  three  survey 
stations,  and  that  the  values  of  wall  temperature  for  the  uncooled 
tests  are  significantly  lower  than  the  recovery  temperature,  TR  , 
for  a  turbulent  boundary  layer  at  Mach  3.  This  characteristic  is  due 
to  heat  losses  from  the  tunnel  to  the  surrounding  room  and,  as  a  result, 
the  data  obtained  with  zero  coolant  flow  had  significant  heat-transfer 
rates . 

typical  velocity  and  total  temperature  profiles  are  shown  in  Fig.  8. 
From  these  profiles  it  can  be  seen  that  wall  cooling  had  little  or  no 
effect  on  the  shape  of  the  flat-plate  velocity  profiles  which  were 
approximately  a  l/6-power  profile.  The  total  temperature  profiles  are 
similar  at  the  outer  edge  of  the  boundary  layer  and  only  near  the  wall 
do  they  show  the  effect  of  cooling.  The  theoretical  total  temperature 
profile  for  the  maximum  cooling  condition  is  also  shown  in  the  lower 
figure.  The  theoretical  value  was  determined  from  the  experimental 
velocity  profiles  making  use  of  the  Crocco  energy  integral  which  form¬ 
ulates  a  linear  relationship  between  total  temperature  and  velocity. 

The  data  for  the  maximum  cooling  condition  exhibit  higher  total  tempera¬ 
tures  than  the  theory.  This  discrepancy  probably  occurs,  in  part, 
because  the  measured  profile  did  not  develop  over  an  isothermal  surface. 

The  initial  2k  in.  of  flat-plate  surface  was  not  cooled  (Tw/TTo«0.85<Tr/TTcj), 
and  a  sudden  drop  in  the  wall  temperature  occurred  at  the  start  of  the 
cooled  flat  plate  (7  in.  upstream  of  the  first  survey  station).  As  a 
consequence,  the  temperatures  in  a  portion  of  the  boundary  layer  away  from 
the  wall  would  be  expected  to  be  higher  than  the  values  for  a  flat-plate 
profile  which  developed  over  a  uniformly  cooled  surface.  The  theoretical 
curve  corresponding  to  the  data  obtained  with  no  coolant  flow  is  not 
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shown  In  Pig.  8  as  it  would  fail  along  the  line  faired  through  the 
uncooled  data  points. 

Experimentally  determined  variations  of  momentum  thickness, 
displacement  thickness,  and  total  enthalpy  thickness  are  shown  in 
Fig.  9  for  each  of  the  three  cooling  conditions.  The  effect  of  wall 
cooling  on  the  integral  parameters  can  be  seen  for  8*  and  6  j  however 
there  is  no  consistent  variation  of  the  boundary  layer  momentum  thick¬ 
ness  with  cooling.  At  the  upstream  station,  the  data  exhibit  an 
increase  in  9  with  increased  cooling,  at  the  middle  station  the  data 
show  little  effect  of  cooling,  and  at  the  downstream  station  the  data 
show  an  opposite  trend.  The  accuracy  in  the  measurements  of  B  is 
estimated  to  be  on  the  order  of  3  to  5#.  Hence  no  significant  effect 

of  cooling  on  momentum  thickness  can  be  deduced  from  the  data  presented 
in  Fig.  9.  presented 


The  theoretical  variations  of  9  ,  8*  ,  and  are  shown  on  Fig.  9 
for  two  values  of  the  wall  temperature  ratio.  The  method  of  calculation 
is  discussed  under  ANALYTICAL  INVESTIGATION.  The  values  of  momentum 
thickness  (0.018  in.)  and  total  enthalpy  thickness  (O.OOh  in.)  at  the 
start  of  the  cooled  portion  of  the  wall  were  estimated  by  extrapolating 
the  data  to  the  start  of  the  cooled  surface  making  use  of  the  theoretical 
slope  to  obtain  a  best  fit.  For  a  wall  temperature  ratio  of  0.8>  (uncooled) 
the  theoretical  and  experimental  variations  of  momentum  thickness  with  x 
are  in  agreement.  The  theoretical  increase  of  the  momentum  thickness  due 
to  cooling  at  the  last  survey  station  is  about  956. 

The  theoretical  values  for  the  displacement  thicknesses  were 
obtained  from  the  theoretical  momentum  thickness  and  the  theoretical 
shape  parameter,  hi/7  ,  for  a  l/7-power  profile  as  tabulated  in  Ref.  33. 

The  displacement  thickness  data  are  in  good  agreement  with  theory  at  a 
wall  temperature  ratio  of  O.85;  however,  the  decrease  in  displacement 
thickness  with  wall  cooling  is  not  as  great  as  that  predicted  by  Ref  33 
As  shown  in. Fig.  8,  the  experimental  total  temperatures  for  the  maximum^ 
cooling  condition  are  greater  than  the  theoretical  values  obtained  by 


17 


ASD  TDR  62-67 


employing  the  Croeco  energy  Integral.  The  total  temperature*  used 
in  evaluating  the  parameters  tabulated  in  Ref.  33  were  coaputed  by 
a  modification  of  the  Crocco  energy  integral  vhich  results  in  even 
lover  total  temperatures  than  would  be  computed  from  a  linear  total 
temperature-velocity  relationship.  Ibis  difference  between  the 
theoretical  and  experimental  values  of  total  temperature  is  the  primary 
reason  for  the  large  discrepancy  between  the  theoretical  and  measured 
values  of  8*  for  the  intermediate  and  maximum  cooling  conditions . 

The  data  for  total  enthalpy  thickness  shov  the  same  trends  as  the 
theory;  however,  the  v&lues  are  not  in  quantitative  agreement.  One 
reason  for  this  lack  of  agreement  can  be  attributed  tc  the  accuracy 
of  the  4>  i_  '  egration .  Small  errors  in  the  measurement  of  total 
temperature  near  the  edge  of  the  boundary  layer,  where  the  mass  flow 
is  large,  can  result  in  large  errors  in  the  calculated  values  of  total 
enthalpy  thickness. 

Figure  10  shows  the  experimental  values  of  the  shape  parameter,  H, 
and  the  reduced  shape  parameter,  6,  for  the  flat  plate  at  Mach  3.  The 
theoretical  values  of  the  shape  parameter  for  a  boundary  layer  with 
constant  total  temperature  (  h{,7  ,  see  Ref.  33)  and  for  a  boundary 
layer  with  a  total  temperature  profile  determined  by  a  modification 
of  the  Crocco  integral  (  H)/7  ,  see  Ref.  33)  are  also  shown  in  thiB 
figure.  Although  the  experimental  shape  parameter  is  reduced  with 
cooling,  the  reduction  is  not  as  great  as  that  predicted  by  Ref.  33. 

This  departure  from  theory  can  also  be  seen  in  the  measured  displace¬ 
ment  thicknesses  discussed  previously. 

Shock  Generator 

The  shock-generator  model  (see  Fig.  l)  was  employed  to  obtain 
boundary  layer  profile  data  immediately  downstream  of  the  shock  reflection 
on  a  cooled  flat  plate  for  generator  angles  of  attack  of  0,  k,  8,  and 
10  deg.  Figure  11  contains  a  schlieren  photograph  and  a  sketch  of  the 
Bhock-ref lection  interaction  region  showing  the  incident  shock  wave, 
the  reflected  shock  wave  (an  initial  compression  followed  by  an 
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expansion  and  then  recompression),  and  the  expansion  fan  originating 
at  the  base  of  the  shock  generator  wedge.  The  inviseid  shock-reflec¬ 
tion  point,  xs  ,  and  the  probe  station  are  also  shown  on  this  figure. 
The  location  of  the  surface  of  the  flat  plate  in  the  schlieren  photo¬ 
graph  was  obscured  by  scratches  in  the  tunnel  window  and  is  therefore 
indicated  by  the  white  line. 

The  measured  wall  static  pressure  and  wall  temperature  distribu¬ 
tions  are  given  in  Fig.  12  for  each  of  the  generator  angles  for  both 
zero  and  maximum  coolant  flow  rates,  hereafter  referred  to  as  "uncooled" 
and  "cooled".  The  location  of  the  inviscid  shock-reflection  point 
varied  from  xs  *  3.3  in.  for  the  10-deg  shock  to  xs  *  3.6  in.  for  the 
O-deg  wave.  The  probe  station  was  at  x  *  4.0  in.  for  all  profiles. 

For  all  generator  angles  of  attack,  the  interaction  started  farther 
downstream  with  wall  cooling  then  without,  and  the  static  pressure 
rose  to  a  higher  value  at  the  probe  station.  This  is  in  agreement 
with  the  trends  shown  by  the  data  of  Ref.  3.  At  a  generator  setting 
of  0  deg,  a  weak  wave  was  observed  in  the  schlieren  even  though  the 
net  flow  deflection  was  very  small.  A  slight  rise  in  wall  static 
pressure  can  be  seen  at  this  generator  setting.  The  pressure  ratios 
at  the  probe  station  were  less  than  the  theoretical  values  for  angles 
of  8  and  10  deg  because  the  expansion  fan  from  the  wedge  intersected 
the  boundary  layer  before  the  pressure  rise  was  completed. 

Figures  13  and  14  present  the  velocity  and  total  temperature 
profiles  at  the  probe  station  for  tests  with  both  cooled  and  uncooled 
walls.  In  Fig.  13  the  velocities  are  ratloed  to  the  nominal  tunnel 
velocity,  and  therefore  do  not  necessarily  approach  unity  at  the  edge 
of  the  boundary  layer.  At  zero  angle  of  attack  the  velocity  profiles 
were  similar  (no  significant  change  in  the  velocity  profile  would  be 
expected  since  the  flat-plate  profiles  showed  no  significant  change 
with  cooling);  at  4  deg  the  velocity  profile  was  much  fuller  for  the 
maximum  cooling  condition;  at  8  deg  the  profiles  were  similar;  and  at 
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10  deg  the  velocity  profile  wi  fuller  for  the  uhbooled  testa.  This 
reversal  of  cooling  effect  with  increasing  angle  of  attack  can  he 
attributed  to  two  opposing  factors:  (l)  for  a  given  pressure  rise 
cooling  tends  to  laprove  the  boundary  layer  profile  since  upstreaa  of  the 
interaction  the  cooled  profile  has  higher  Mach  numbers  through  the 
boundary  layer  and  is  able  to  negotiate  the  pressure  rise  with  less 
reduction  in  velocity  and  hence  less  distortion  than  the  uncooled  pro¬ 
file,  and  (2)  for  a  given  incident  shock  strength  cooling  tends  to 
(a)  reduce  the  length  of  the  interaction  (i.e.,  larger  adverse  pressure 
gradient)  and  (b)  shift  the  interaction  downstream,  effectively  shifting 
the  probe  station  forward  into  the  highly  distorted  portion  of  the 
interaction  region.  These  factors  tend  to  Increase  the  distortion 
measured  at  the  probe  station.  From  the  data,  it  appears  that  the  first 
of  these  factors  predominates  at  low  angles  of  attack,  while  at  higher 
angles  of  attack,  the  second  of  these  factors  was  more  pronounced. 

The  boundary  layer  momentum  and  displacement  thicknesses  and  the 
reduced  shape  parameters  for  these  profiles  are  presented  in  Fig.  15. 

The  figure  shows  the  expected  trend  with  increasing  angle  of  attack, 
l.e.,  the  values  of  Q ,  8*  ,  and  G  generally  Increase  with  Increasing 
shock  strength.  The  Increase  in  G  between  0-  and  10-deg  angles  of 
attack  is  greater  for  the  data  with  maximum  coolant  flow  than  for  the 
data  with  zero  coolant  flow.  This  is  due  to  the  same  factors  which 
produced  the  reversal 'in  trend  in  the  velocity  profiles.  No  consistent 
trend  could  be  established  for  the  corrected  total  enthalpy  thickness 
data  obtained  with  the  shock  generator  and,  consequently,  these  data 
are  not  presented. 

Circular  _Arc_ 

Boundary  layer  profile  data  were  obtained  at  five  stations  on 
the  circular-arc  model  for  three  rates  of  coolant  flow  which  corres¬ 
ponded  to  wall  temperature  ratios  of  approximately  0.85,  0.60,  and 
0.45  (the  uncooled,  intermediate  cooling,  and  maximum  cooling  conditions, 
respectively).  Schlieren  photographs  of  the  flow  over  the  circular-arc 
model  with  each  of  these  cooling  rates  are  shown  in  Fig.  16.  Figure  17 
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staova  the  measured  and  theoretical  values  of  nil  static  pressure,  , 
the  calculated  static  pressure  at  the  edge  of  the  boundary  layer,  W, 
and  the  measured  wall  temperature,  Tw  ,  as  a  function  of  distance 
the  surface.  The  stations  at  which  surveys  were  obtained  are  noted  on 
the  figure.  The  vail  teaperature  ratio  did  not  vary  significantly  over 
the  portion  of  the  aodel  covered  by  the  survey  stations,  as  shown  by  the 
lover  portion  of  the  figure. 

The  velocity  profiles  obtained  on  the  circular-arc  aodel  are  shown 
in  Pig.  18.  The  velocity  profiles  upstream  of  the  tangent  point  show 
no  significant  change  with  cooling,  as  would  be  expected  from  the  flat- 
plate  data,  whereas  the  velocity  profiles  on  the  curved  portion  of  the 
surface  do  show  a  substantial  improvement  in  shape  and  a  decrease  in 
boundary  layer  height  with  cooling.  This  trend  with  Increased  coolant 
flow  Is  observed  at  all  stations  except  at  the  12-dep  station  where  the 
velocity  profile  vith  maximum  cooling  shows  negligible  change  from  the 
profile  obtained  at  the  intermediate  cooling  rate.  Analysis  of  the 
data  at  this  station  indicates  that  there  was  an  error  In  the  probe 
position  data  due  to  backlash  in  the  probe  drive  mechanism.  An  attempt 
has  been  made  to  correct  the  data  but  the  magnitude  of  the  correction 
(Ay  *  0.015  in.)  could  not  be  accurately  determined.  The  total- 
temperature  profiles  (Pig.  19)  on  the  upstream  portion  of  the  curved 
surface  show  a  change  with  cooling  over  the  lover  third  of  the  profile 
as  was  the  case  with  the  flat-plate  profiles  (see  Pig.  8).  On  the 
downstream  portion  of  the  surface,  however,  the  profiles  show  no 
change  with  cooling  except  in  the  immediate  vicinity  of  the  wall.  There 
is  a  strong  possibility  that  the  profiles  at  the  downstream  station  could 
be  seriously  affected  by  probe  conduction  error  due  to  the  low  bleed 
flow  through  the  probe  in  the  distorted  profiles  at  these  stations.  For 
a  given  height  within  the  boundary  layer  at  the  24-deg  station,  boundary 
layer  mass  flow  calculations  indicate  that  there  is  approximately  3056 
more  mass  flow  for  the  maximum  cooling  rate  than  there  is  for  the 
uncooled  profile.  This  is  an  important  consideration  in  the  diffusion 
of  a  supersonic  boundary  layer  (e.g.,  such  as  encountered  in  high- 
performance  internal-compression  inlets)  since  this  increase  in  mass 
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flow  In  the  boundary  layer  will  permit  the  use  of  higher  inlet 
contraction  ratios  with  cooled  walls  than  without. 

Figure  20  shows  the  experimental  variations  of  displacement 
thickness  reduced  shape  parameter.  The  change  of  displacement 
thickness  along  the  surface  was  less  with  wall  cooling  than  without; 
and  the  data  for  reduced  shape  parameter  show  that  cooling  produces 
a  significant  reduction  (i.e.,  improvement)  in  the  shape  parameter. 

Figure  21  shows  the  variation  of  boundary  layer  momentum  thickness  and 
corrected  total  enthalpy  thickness  computed  from  the  measured  profiles. 
The  corrected  total  enthalpy  thickness,  <f>c  ,  is  defined  as  the  local  <f> 
times  the  ratio  of  the  mass  flow  computed  from  the  local  wall  static 
pressure  ratio  to  the  mass  flow  at  the  nominal  free-stream  conditions. 
This  parameter  is  more  useful  than  4>  since  it  is  directly  proportional 
to  the  total  heat  removed  from  the  boundary  layer,  i.e.,  Q  =  />o uo ^To • 
The  rate  of  increase  of  momentum  thickness,  9  ,  (top  half  of  Fig.  21) 
over  the  upstream  portion  of  the  surface  is  the  same  for  the  three 
cooling  conditions.  Over  the  downstream  portion  of  the  surface,  the 
rate  of  increase  of  9  is  significantly  reduced  with  wall  cooling  as  would 
be  expected  on  the  basis  of  the  improved  velocity  profiles  shown  in 
Fig.  10.  The  corrected  total  enthalpy  thickness  increases  with  cooling 
at  each  station,  but  there  is  a  lack  of  consistency  in  the  variation  with 
distance  over  the  downstream  portion  of  the  plate.  This  may  also  be 
attributed  to  probe  conduction  error  since  such  an  error  would  tend  to 
increase  <£c  for  the  case  without  cooling  and  reduce  <j>c  for  the  case 

with  cooling. 

Theoretical  values  for  displacement  thickness,  reduced  shape 
parameter,  momentum  thickness,  and  corrected  total  enthalpy  thickness 
are  also  shown  with  the  experimental  data  in  Figs.  20  and  21.  The 
theoretical  values  are  based  on  the  methods  discussed  in  a  previous 
section  (ANALYTICAL  INVESTIGATION ) .  For  each  value  of  wall  temperature 
ratio,  the  initial  value  of  G  on  the  circular  arc  was  determined  from  the 
Mach  3  flat-plate  studies  {see  Fig.  10).  The  theoretical  values  of  G  (see 
Fig.  20)  are  in  general  agreement  with  the  experimental  data  except  in 
the  uncooled  case  where  the  theoretical  values  are  a  maximum  of  10$  lower 
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STS5  SS  S  the 'experimental  data.  The  largest  difference. 
Fig.  are  .  for  ^  uncooled  case  where  the 

between  theory  and  “P«‘  Btltlcn  are  30  and  35*  of  the  experl- 
differences  at  the  last  Surrey  s  agreement  between 

:^y  bT^cte/for  the  uncooled  case  since  the 

lag- length  theory  expired  Thef.  1 
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most  of  the  data  fell  with!  35j  (  differences  shown  in  Figs.  20 
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calculation  procedure.  As  can  be  M«  **£■  "ter, 

between  the  theoretical  and  experimental  values  or  tne  s  ,  par 
c  ,  are  of  the  order  of  the  experimental  accuracy  ( »  3* >  at  £1  ***“ 
’  t  the  two  points  obtained  without  cooling  at  the  12- and  10  <*g 
Pt  _  P  .  d  dflta  were  expected  to  show  the  best  agreement 
^^rv^nce  th^Sry  was  originally  developed  for 
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on  the  method  of  analys  .  lac- length  procedure  and 
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Of  the  lag- length  procedure  alone  (see  Ref.  lj. 

ae  lower  portion  of  Fig.  21  shows  the  theoretical  values  for  the 
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of  the  oints  at  the  M,,  and  possibly  18-dcg  statlons  cou  ^ 
seriously  affected  by  probe  conduction  error.  A  statistical  a. caly 
of  the  data  indicated  that  the  value  of  total  heat  transfer  is  greater 
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than  the  lower  probable  limit  by  an  amount  which  is  approximately 
30^  of  the  difference  between  the  upper  and  lower  probable  limits. 

Mach  6  Studies 

Flat_Plate__ 

Tests  were  conducted  at  Mach  6  with  the  flat-plate  model  installed. 
The  wall  Mach  number  and  wall  temperature  distributions  for  zero  coolant 
flow  rate  (uncooled)  and  maximum  coolant  flow  rate  (cooled)  are  snown 
in  Fig.  22.  The  Mach  number  varied  from  a  value  of  6.02  at  the  up¬ 
stream  survey  station  to  a  value  of  5*®5  at  the  downstream  station; 
the  wall  temperature  was  uniform  over  the  three  survey  stations.  As 
with  the  Mach  3  data,  the  wall  temperature  for  the  uncooled  tests  was 
lower  than  the  turbulent  recovery  temperature;  with  maximum  cooling, 
however,  the  wall  temperatures  were  lower  than  for  the  Mach  3  tests. 
Typical  velocity  and  total  temperature  profiles  for  the  two  cooling  rates 
are  shown  in  Fig.  23.  The  velocity  profiles  were  similar  (approximately 
a  l/5-power  profile);  the  cooled  profile  was  a  little  fuller  near  the 
wall.  As  in  the  Mach  3  studies,  the  total  temperatures  near  the  edge 
of  the  boundary  layer  were  the  same  and  the  effect  of  cooling  was 
evident  in  the  lower  portion  of  the  boundary  layer.  The  theoretical 
total  temperature  distribution  for  a  linear  temperature-velocity 
relationship  as  formulated  by  the  Crocco  energy  integral  is  shown  for 
the  profile  obtained  with  the  maximum  coolant  flow  rate.  Again,  the 
experimental  total  temperatures  are  significantly  different  from  the 
theoretical  values  (see  Mach  3  flat-plate  results);  a  characteristic 
which  may  be  attributed  to  the  stepwise  temperature  distribution  on 
the  surface  over  which  the  boundary  layer  developed.  As  in  the  case 
of  the  Mach  3.0  flat  plate  data,  the  theoretical  total  temperature 
profile  with  no  cooling  agrees  with  the  total  temperature  data  points. 

The  theoretical  and  experimental  total  temperature  distributions  are 
in  agreement  for  the  cases  with  no  cooling  because  there  is  no  step¬ 
wise  temperature  gradient  along  the  surface. 

The  experimentally  obtained  variations  of  boundary  layer  momentum 
thickness,  displacement  thickness,  and  total  enthalpy  thickness  are 
shown  in  Fig.  24  together  with  the  theoretical  values  for  these  para- 
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ae ter  1 .  lhe  data  show  a  larger  momentum  thickness,  a  smaller  dis¬ 
placement  thickness,  and  a  larger  total  enthalpy  thickness  for  the 
cooled  profiles  for  the  uncooled  profiles.  The  direction  of  all 
these  trends  is  In  general  agreement  with  theory.  The  values  of 
momentum  thickness  (0.030  in.)  and  total  enthalpy  thickness  (0.010  in.) 
at  the  start  of  the  cooled  portion  of  the  surface  were  estimated  by 
extrapolating  the  data  to  the  start  of  the  cooled  surface.  The  momentum 
thickness  data  are  in  general  agreement  with  the  theoretical  values; 
the  deviation  is  of  the  order  of  the  repeatability  of  the 

measurements  as  Indicated  by  the  two  points  for  the  uncooled  tests  at 
the  middle  survey  station,  lhe  change  in  momentum  thickness  from  the 
uncooled  case  to  the  cooled  case  is  in  good  agreement  with  theory  but 
is  of  the  same  order  of  magnitude  as  the  accuracy  of  the  data;  hence  no 
significant  conclusion  may  be  drawn  from  this.  The  data  show  higher 
displacement  thicknesses  than  predicted  by  theory  and,  similar  to  the 
Mach  3  data,  the  reduction  in  displacement  thickness  with  cooling  is  not 
as  great  as  that  predicted  by  theory.  The  data  for  total  enthalpy 
thickness  are  in  fair  agreement  with  the  theoretical  values  with 
mviimiin  cooling,  but  the  data  obtained  on  the  uncooled  surface  do  not 
agree  with  theory.  As  was  pointed  out  for  the  Mach  3  studies,  the 
probable  source  of  error  is  in  the  measurement  of  total  temperature 
near  the  edge  of  the  boundary  layer  which  can  result  in  large  errors  in 
total  enthalpy  thickness  due  to  the  high  mass  flow  rates  at  the  edge  of 
the  boundary  layer. 

The  variation  of  the  boundary  layer  shape  parameters  with  wall 
temperature  ratio  is  shown  in  Fig.  25.  The  theoretical  values  of  the 
shape  parameters  for  cooled  1/5-  and  l/7-power  profiles  were  obtained 
from  Ref.  33.  The  experimental  values  of  the  shape  parameters  are 
reduced  by  wall  cooling,  but  the  reduction  is  not  as  great  as  that 
predicted  by  the  temperature -velocity  relationship  for  the  modified 
Crocco  integral  as  tabulated  in  Ref.  33- 

I s  entropic_Surface_ 


Boundary  layer  surveys  were  obtained  at  five  stations  on  the 
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ltentropic-surfkce  model  for  zero  and  maxima  coolant  flow  rates.  The 
vail  temperature  ratios  were  approximately  0.82  (uncooled)  and  0.1*5 
(cooled).  Schlleren  photographs  of  the  flow  over  the  lsentroplc-surffcee 
model  at  each  of  these  temperature  ratios  are  shown  in  Fig.  26.  Measured 
and  theoretical  distributions  of  vail  static  pressure,  measured  vail 
temperature,  and  the  calculated  distribution  of  static  pressure  at  the 
edge  of  the  boundary  layer  are  shown  in  Fig.  27 •  Good  agreement  was 
obtained  between  the  measured  and  theoretical  values  of  vail  static 
pressure.  There  was  no  significant  change  in  vail  temperature  along 
the  surface  of  the  model  for  the  data  obtained  vith  the  uncooled  vail, 
but  the  data  for  the  cooled  wall  show  a  fairly  large  longitudinal  tem¬ 
perature  gradient.  Die  data  for  this  case  are  reported  at  an  average 
temperature  ratio  of  0.1*5  (O.36  at  the  upstream  survey  station  and  0.57 
at  the  last  survey  station).  The  survey  stations  and  values  of  local 
surface  angles  are  shown  in  the  lover  portion  of  the  figure.  Although 
the  model  had  a  total  turning  angle  of  3I*  deg,  the  total  turning  *nei#» 
at  the  last  station  at  which  surveys  were  obtained  was  26.5  deg.  Down¬ 
stream  of  this  station  the  flow  was  influenced  by  the  expansion  orig¬ 
inating  at  the  end  of  the  curved  surface.  On  this  portion  of  the 
surface,  the  turning  angle  increases  from  26.5  to  34.0  deg  in  a  dis¬ 
tance  of  0.31  in.  whereas  the  boundary  layer  thickness  was  0.4  in. 

Velocity  and  total  temperature  profiles  are  shown  in  Figs.  28  and 
29.  The  velocity  profiles  at  the  1-deg  station  show  no  significant 
change  with  cooling,  whereas  the  velocity  profiles  at  the  downstream 
stations  show  a  substantial  improvement  in  shape  and  a  decrease  in 
boundary  layer  height  with  cooling.  The  total  temperature  profiles 
at  the  1-deg  station  appear  to  be  much  the  same  as  the  flat-plate  total 
temperature  profiles  with  almost  all  of  the  change  due  to  cooling 
occurring  over  the  lower  quarter  of  the  profile.  Downstream,  the 
profiles  differ  farther  away  from  the  wall  arid,  at  the  12.8-deg  station, 
the  difference  is  seen  over  half  of  the  profile.  For  the  last  two 
stations  on  the  surface  the  measured  total  temperature  profiles  for 
the  cooled  case  are  higher  than  those  for  the  uncooled  case  in  the 
middle  third  of  the  profile.  In  spite  of  this  apparent  discrepancy, 
the  total  temperature  for  the  cooled  case  is  always  less  than  or  equal 
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to  the  total  temperature  for  the  uncooled  case  when  total  temperature 
le  plotted  as  a  function  of  the  Integrated  mass  flow  since  the  cooled 
profile-3  have  more  mast  flow  near  the  wall  than  the  uncooled  profiles. 
Calculations  of  the  mass  flow  in  the  boundary  layer  show  that  the 
profiles  at  the  24-deg  station  for  the  cooled  case  have  500  more  mass 
flow  than  the  uncooled  profiles .  As  was  the  case  for  Mach  3,  this 
increased  mass  flow  within  a  given  height  implies  that  supersonic 
diffusers  having  cooled  walls  could  operate  with  higher  area  contrac¬ 
tion  ratios  and  hence  higher  total  pressure  recoveries  than  uncooled 
inlets . 

The  data  for  displacement  thickness  and  reduced  shape  parameter 
are  presented  in  Fig.  30.  The  large  decrease  in  8*  for  both  cases 
is  due  to  the  large  compression  of  the  boundary  layer  over  the  length 
of  the  surface.  Cooling  caused  an  additional  decrease  in  the  dis¬ 
placement  thickness  at  the  downstream  stations.  At  all  stations 
cooling  caused  a  significant  reduction  (improvement)  in  the  shape 
parameter.  Figure  31  shows  the  variation  of  momentum  thickness  and 
corrected  total  enthalpy  thickness  with  distance  along  the  surface. 

The  difference  in  initial  level  of  the  two  curves  is  of  the  same  order 
of  magnitude  as  the  accuracy  of  the  data  and  hence  is  not  significant; 
however,  the  over-all  reduction  in  9  for  the  cooled  surface  is  greater 
than  that  for  the  uncooled  case  as  would  be  expected  on  the  basis  of 
the  improved  velocity  profiles  shown  in  Fig.  28.  The  corrected  total 
enthalpy  thickness  data  show  a  consistent  increase  with  wall  cooling. 

The  theoretical  predictions  of  displacement  thickness,  reduced 
shape  parameter,  momentum  thickness,  and  the  upper  and  lower  probable 
limits  of  corrected  total  enthalpy  thickness,  are  also  shown  in  Figs. 

30  and  31.  These  predictions  are  based  on  the  methods  outlined  in  a* 
previous  section  (see  ANALYTICAL  INVESTIGATION).  The  theoretical  values 
of  the  reduced  shape  parameter  (Fig.  30)  show  the  same  trends  as  the 
data,  but  the  quantitative  agreement  is  not  as  good  as  it  was  for 
Mach  3  circular-arc  studies .  The  initial  values  of  G  were  determined 
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from  the  Mach  6  flat-plate  data  (see  Fig.  25).  As  can  be  seen  from 
Fig.  30  the  average  deviation  between  the  theoretical  and  experimental 
values  of  the  shape  parameter  is  on  the  order  of  25$.  As  shown  in 
APPENDIX  III,  the  errors  for  the  Mach  6  data  introduced  by  the  use  of 
Eq.  8  are  less  than  100,  thus  indicating  that  the  errors  introduced 
by  the  lag- length  procedure  are  on  the  order  of  20  to  250  for  the  Mach  6 
data.  This  Indicates  the  need  for  additional  data  for  use  in  extending 
the  lag- length  correlation  to  higher  Mach  numbers .  The  data  on  which 
the  existing  correlation  is  based  (see  Ref.  l)  were  obtained  in  the 
Mach  number  range  of  2.5  to  3*5 •  The  theoretical  values  of  the  dis¬ 
placement  thickness  (top  of  Fig.  30)  in  agreement  with  the  data; 
however,  the  theoretical  values  of  momentum  thickness  over  the  down¬ 
stream  portion  of  the  surface  (see  Fig.  31)  significantly  higher 
than  the  data.  The  maximum  difference  is  about  350  of  the  measured 
value.  This  difference  is  the  same  order  of  magnitude  as  that  obtained 
for  the  Mach  3  circular-arc  studies;  but  it  is  in  the  opposite  direction. 
For  both  configurations  (Mach  3  circular  arc  and  Mach  6  lsentropic  sur¬ 
face),  the  maximum  discrepancies  are  within  the  limits  of  accuracy  for 
the  lag-length  calculation  procedure  developed  in  Ref.  1.  It  should  be 
noted,  however,  that  the  data  on  which  the  lag-length  theory  is  based 
were  obtained  in  the  Mach  number  range  of  2.0  to  3.5.  An  alternate 
calculation  of  the  theoretical  boundary  layer  growth  over  the  Mach  6 
isentropic  compression  surface  was  made  by  employing  the  method  devel¬ 
oped  in  Ref.  5.  This  calculation  predicted  boundary  layer  separation 
after  18  deg  of  isentropic  compression,  whereas  the  experimental 
results  showed  that  the  boundary  layer  did  not  separate  at  any  point 
on  the  isentropic-corapression  surface  having  a  total  turn  of  34  deg. 

Thus,  even  though  there  are  relatively  large  differences  between  the 
lag- length  theory  and  experiment,  the  values  obtained  from  the  lag- 
length  theory  are  in  better  agreement  with  experiment  than  other  avail¬ 
able  theories  for  predicting  supersonic  turbulent  boundary  layer  growth 
over  cooled  walls  in  adverse  pressure  gradients. 

The  theoretical  predictions  for  the  upper  and  lower  probable  limits 
of  corrected  enthalpy  thickness  are  presented  in  the  lower  portion  of 
Fig.  31.  For  the- cooled  case,  the  data  agree  very  well  with  the  values 
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for  the  lover  limit  down  to  the  5. 5 -deg  station  and  for  downstream 
stations  tend  to  fall  between  the  two  limits.  For  the  uncooled  case, 
the  data  also  fall  along  the  lover  probable  limit  to  the  5-5-deg 
station,  but  are  above  the  upper  limit  downstream  of  the  12.8-deg 
station.  Here  again,  there  is  a  possibility  of  probe  conduction 
error  u.it  to  the  low  bleed  flow  through  the  probe  in  the  distorted 
profiles  at  the  two  downstream  stations.  Complete  elimination  of 
this  error  would  probably  tend  to  reduce  the  in  the  uncooled 
case  and  increase  4>c  in  the  cooled  case.  A  statistical  analysis  of 
the  data  has  indicated  that  the  total  heat  flow  is  greater  than  the 
lover  probable  limit  by  an  amount  which  is  approximately  30$  of  the 
difference  between  the  upper  and  lover  probable  limits. 


CORRELATION  OF  SKIN-FRICTION  AND 
HEAT-TRANSFER  COEFFICIENTS 


Calculation  of  heat  transfer  in  both  subsonic  and  supersonic 
flows  is  in  general  based  on  a  form  of  Reynolds  analogy  which  relates 
the  Stanton  number  (heat-transfer  coefficient),  St  ,  to  the  skin- 
friction  coefficient,  Cf  •  This  relationship  can  be  evaluated  with 
the  aid  of  the  Crocco  energy  integral  (see  Refs.  3^  an d  36).  In  its 
simplest  form  the  Crocco  energy  integral  can  be  written  as  follows: 


CpTj  s  Cj  + CgU 


(ID 


in  which  C(  and  C2  are  evaluated  from  the  boundary  conditions  at  the 
wall  where  Tt=  Tw  and  u  =0  and  in  the  free  stream  where  Tt  =  TTq  and 
u  =  u0 .  It  should  be  noted  that  the  Crocco  integral  was  derived  for 
laminar  flow  and  for  a  Prandtl  number  of  unity;  hence  its  application 
to  the  turbulent  flow  of  a  gas  having  a  nonunity  Prandtl  number  is 
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questionable.  However,  the  type  of  data  presented  herein  (velocity 
and  total- temperature  profiles)  allows  a  direct  check  on  the  applica¬ 
bility  of  the  Crocco  integral  to  flows  of  the  type  studied.  Data 
from  one  of  the  Mach  6,  cooled,  flat-plate  profiles  are  cross  plotted 
in  terms  of  total  temperature  and  velocity  in  Pig.  32*  la  this 
representation  the  value  for  the  Crocco  integral  would  be  a  straight 
line  between  the  wall  temperature  ratio  (at  aero  velocity)  and  unity 
total  temperature  ratio  (at  unity  velocity  ratio).  It  can  be  seen 
from  this  figure  that  the  measured  total  temperature  deviates  from 
the  values  for  the  Crocco  integral  especially  for  the  higher  velocity 
ratios  where  it  is  higher  than  that  predicted  by  the  Crocco  Integral. 
These  higher  total  temperatures  may  be  a  result  of  the  stepwise  tempera¬ 
ture  distribution  (uncooled  over  initial  portion,  cooled  over  rear 
portion)  on  the  flat  plate  over  which  the  boundary  layer  developed. 

Most  of  the  mass  flow  in  the  boundary  layer  is  in  the  high  velocity 
region;  hence,  the  momentum  thickness,  Q  ,  and  the  total  enthalpy 
thickness,  <f>  ,  obtained  by  integrating  the  profile  data  are  appre¬ 
ciably  different  from  those  predicted  on  the  basis  of  a  Crocco  energy 
integral . 

All  the  flat-plate  velocity  profiles  were  plotted  on  temperature- 
velocity  coordinates  in  the  manner  of  Pig.  32  and,  in  general,  a 
straight  line  drawn  between  the  wall  temperature  and  the  data  points 
in  the  lower  velocity  region  of  the  boundary  layer  has  the  same  slope 
as  the  data  points  in  that  region.  It  should  be  noted  that  the  first 
data  point  having  a  non-zero  velocity  was  obtained  with  the  probe 
touching  the  wall  and  hence  the  measured  temperature  was  in  error  as 
a  result  of  heat  conduction  from  the  probe  to  the  wall.  Consequently, 
no  attempt  was  made  to  fair  the  line  through  the  first  data  point  in 
the  profile.  The  linear  total- temperature  relationship  over  the  lower 
velocity,  portion  of  the  profile  suggested  that  locally  (in  the  laminar 
sublayer)  a  Crocco  energy  integral  is  valid  and  hence  the  temperature- 
velocity  derivative  at  the  wall  is  adequately  defined  and  can  be 
employed  in  evaluating  the  ratio  of  the  Stanton  number  to  skin-friction 
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coefficient,  St/Cf  .  It  should  be  noted  that  a  direct  evaluation  of 
the  heat  transfer  rate  (k  &  |  w  )  and  wall  shearing  stress  ( /*  &U- 
could  not  be  Bade  from  the  derivatives  of  the  temperature  and  velx>£ jft;  rw 
with  respect  to  y  because  a  straight-line  fairing  frost  the 
condition  did  not  blend  with  the  profile  data  near  the  wall. 


) 


The  following  equation  was  employed  in  relating  the  Stanton 
to  the  skin-friction  coefficient  (see  APPENDIX  II  for  derivation  of 
derivative  method). 


(12) 


In  order  to  evaluate  Eq.  12  the  tenrperature-velocity  derivatives  at  the 
wall  were  plotted  as  a  function  of  the  difference  between  the  recovery 
and  measured  wall-temperature  ratios,  as  shown  in  Pig.  33.  When  plotted 
in  this  form  the  data  had  a  linear  relationship  with  very  little  scatter 
and  there  was  no  significant  difference  between  the  Mach  3  and  Mach  6 
flat-plate  data.  A  least-squares  line  originating  at  the  origin  was 
faired  through  all  the  data  to 'establish  an  average  slope  for  evaluating 
St/Cf  by  Eq.  12.  For  a  Prandtl  number  of  0.70  (average  value  for  range 
of  wall  temperatures  covered)  the  ratio  of  Stanton  number  to  skin-friction 
coefficients  from  Eq.  12  is  O.95,  with  a  standard  deviation  of  an  indi¬ 
vidual  point  being  about  6 i>.  This  value  (0.95)  is  significantly 
than  that  obtained  from  the  Reynolds -Colburn  analogy  (Ref.  34)  as  follows: 


cj~  =  1p^7T:  0  635  F0R  Pr  5  0  70  (13) 


In  order  to  check  the  order  of  magnitude  of  St/Cf  ,  an  alternate 
method  was  employed.  This  method  used  the  integral  quantities  of 
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■omntuB  thickness,  8  ,  and  total  anthalpy  thickness,  4>  .  It  la 
shown  in  APPENDIX  II  that  tha  followln*  equation  ralates  St  and  Cf 
with  $  and  8  for  flat-plats  boundary  laysr  growth  over  an  isothsnal 
aurfact. 


Si  - 

cf  *  tr  tw  \  (1U) 

'  Tr0  tt0  ' 


Since  the  model  employed  for  the  experiments  did  not  have  an  isothermal 
surface  over  the  entire  length  of  boundary  layer  growth  (the  front  portion, 
about  half  the  total  length,  was  uncooled  and  resulted  in  a  stepwise 
temperature  distribution),  the  effective  wall  temperature  (a  length- 
weighed  average  defined  below)  was  used  in  evaluating  Eq.  14. 


V'  + V*  (15) 

T w  2  — L|  +  i2 

In  this  equation  the  subscript  i  refers  to  the  uncooled  portion  and  the 
subscript  2  refers  to  the  length  of  cooled  portion  up  to-the  measuring 
station.  This  approximation  for  the  effective  temperature  could  be 
improved  by  incorporating  an  average  skin-friction  coefficient  (see 
APPENDIX  II);  however,  the  increased  accuracy  was  not  Justified  on  the 
basis  that  the  total  enthalpy  thickness  integral,  ,  was  already 
believed  to  be  the  least  accurate  quantity  in  the  expression.  The 
ratios  of  <p/8  for  all  the  flat-plate  data  were  plotted  as  a  ^function 
of  the  difference  between  the  recovery  and  effective  wall  temperature 
ratios  as  shown  in  Fig.  34.  A  least-squares  straight  line  was  faired 
through  the  data  and  the  origin.  The  standard  deviation  of  a  point  was 
about  17$,  and  again  there  was  no  significant  difference  between  the 
Mach  3  and  Mach  6  data.  The  ratio  of  St/Cf  evaluated  by  this  integral 


32 


ASD  TER  62-67 


Mtfaod  mi  0.80$.  Although  this  value  of  the  ratio  of  Stanton  number 
to  •kin-friction  coefficient  1*  lower  than  that  obtained  by  the  fir«t 
(differential)  Method,  it,  too,  i»  ■ignificantly  higher  than  the 
theoretical  value  obtained  froa  Eq.  13*  Both  aethode,  however,  are 
•ubject  to  error*  in  the  Measured  total  temperature.  If  the  true 
temperature s  were  lower  than  the  value*  used,  the  ratio  of  St/Cf 
computed  from  the  derivative  Method  would  be  lower  and  the  ratio  of 
St/Cf  computed  by  the  integral  method  would  be  higher  than  the 
valueB  obtained  in  the  present  analysis.  In  order  to  decrease  the 
ratio  of  Stanton  number  to  skin- friction  coefficient  by  both  methods, 
the  true  total  temperature  would  have  to  be  lower  than  the  value  used 
near  the  wall  and  higher  than  the  value  used  near  the  free  stream. 

The  method  used  to  evaluate  total  temperature  (discussed  in  APPENDIX  I) 
has  already  favored  the  data  in  the  direction  of  lower  values  of  St/Cf 
over  more  conventional  ways  of  evaluating  total  temperature  from  the 
temperatures  measured  with  a  shielded  thermocouple. 

A  factor  which  may  contribute  to  the  measurement  of  a  high  value 
for  St/Cf  is  the  use  of  a  wall  which  was  not  cooled  from  the  origin 
of  the  boundary  layer  growth;  heat  transfer  rates  greater  than  the 
flat-plate  value  would  be  expected  at  the  beginning  of  the  cooled 
portion  of  a  plate  having  a  stepwise  temperature  distribution.  This 
problem  is  discussed  in  Ref.  37-  From  the  data  presented  in  this 
reference  the  error  in  heat-transfer  coefficient  that  would  be  expected 
from  the  first  (differential)  method  would  range  from  10  to  15 i,  and 
the  error  expected  from  the  second  (integral)  method  would  range  from 
15  to  25$.  Thus,  errors  due  to  a  stepwise  temperature  distribution 
can  be  quite  large  for  the  present  experiments  and  can  account  for  the 
difference  between  the  theoretical  value,  St/Cf  =  O.635,  and  the 
measured  value,  St/Cf  *  0.8l,  obtained  from  the  integral  method; 
however,  these  errors  do  not  account  for  the  large  difference  between 
the  theory  and  the  value  obtained  by  the  differential  method, 

St/Cf  =  0.95,  which  gave  more  self-consistent  results.  The  value 
calculated  by  the  differential  method  was  based  on  wall  derivatives 
obtained  from  straight-line  fairings  on  the  total  temperature-velocity 
plots  (see  Fig.  32).  More  data  in  the  lower  velocity  region  of  the 


33 


* 


ASD  TUI  62-67 


boundary  layer  is  needed  to  varify  that  such  a  Umc  relationship  la 
valid. 


ginea  tha  derivative  aethod  of  evaluating  tha  patio  of  Stanton 
to  skin-friction  eoaffldant  for  tha  flat-plate  data  gives  the 
■ore  conalatent  results,  thii  aethod  was  alao  alloyed  to  evaluate 
St/Cf  for  the  boundary  layer  growth  In  adverae  preaaure  gradient*. 
The  data  were  evaluated  at  each  aeaauring  atation  on  the  circular-arc 
and  isentropic-surface  models  and  for  each  deflection  angle  of  the 
shock-generator  teata.  The  results  are  presented  In  Fig.  35  *a  a 
function  of  the  total  deflection  angle,  ♦  ,  through  which  flow  vaa 
turned.  The  flagged  symbols  for  the  high  deflection  angles  indicate 
that  St/Cf  may  be  higher  than  that  shown  by  the  data.  For  these 
data  points,  a  straight-line  fairing  on  the  tesrperature-velodty  plot 
(see  Fig.  32)  from  the  wall  did  not  blend  with  the  data  points  in  the 
lover  part  of  the  profile.  The  resulting  error  is  believed  to  be  in 
a  direction  which  reduces  the  measured  values  of  St/Cf  .  The  data 
presented  In  Fig.  35  show  that  the  ratio  of  St/Cf  increases  as  the 
total  deflection  angle  (the  amount  of  deceleration)  is  Increased  and 
that  the  deviation  from  the  flat-plate  value  is  greatest  for  the 
moat  severe  pressure  gradients  (shock  generator)  and  lowest  for  the 
mildest  pressure  gradients  (isentropic  surface).  This  increase  in 
St/Cf  does  not  imply,  however,  that  the  heat  transfer  is  increasing 
at  these  points  since  the  skin-friction  coefficient  is  less  than  the 
flat-plate  value  due  to  the  distorted  profile. 
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APFEKDIX  I 

METHOD  OP  EVALUATING  TOTAL  TEMPERATURE 
THROUGH  A  BOUNDARY  LAYER 


The  temperature  measured  by  a  shielded  thermocouple  probe  does  not 
equal  the  local  total  temperature  of  the  stream.  For  a  shielded  thermo¬ 
couple  probe  of  the  typo  used  herein  (see  Fig.  5)>  the  velocity  of  the 
air  past  the  thermocouple  Junction  is  in  the  very  low  subsonic  speed 
range  and,  as  a  result,  the  recovery  temperature  at  the  Junction  should 
be  within  1/10  of  1 $  of  the  true  total  temperature.  However,  most  probes 
of  this  type  indicate  considerably  lower  temperatures .  The  conventional 
way  of  correcting  the  measured  probe  temperature  is  to  employ  a  recovery 
factor  for  the  probe  based  on  a  probe  calibration.  Most  probes  have 
recovery  factors  varying  from  0.94  to  0.98  (see  Ref.  38)  depending  on 
the  method  of  fabrication.  The  difference  between  the  true  and  measured 
temperatures  is  commonly  cr  ilted  to  radiation  effects  which  has  led  to 
the  use  of  multiple  shields .  However,  after  an  examination  of  the  order 
of  magnitude  of  the  radiation  effect  in  the  probes  employed  herein  it  was 
concluded  that  radiation  was  not  the  pi  imary  cause  of  probe  temperature 
errors ;  the  real  source  of  error  was  felt  to  be  due  to  conduction  from  the 
thermocouple  wires  to  the  average  temperature  of  the  probe  head.  If  the 
error  Is  assumed  to  be  due  to  conduction,  the  following  heat  balance  can 
be  written: 


heat  into  Junction  q  s  St  cp(TT  -  TM)A,  (l6) 

heat  out  of  Junction  q  =  K  j*  (TM  -  Tp)  (17) 

*2 

where  Tt  is  the  local  total  temperature,  TM  the  temperature  of  the 
Junction  (measured  temperature);  Tp  the  average  probe  temperature; 
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A,iAr<t  are  seme  effective  ami  and  lengths  (all  unknown)  which  are 
introduced  to  provide  the  correct  units  j  and  the  bars  (  “  )  denote  average 
conditions.  The  quantities  that  vary  during  a  boundary  layer  survey  are  the 
temperatures  and  the  mss  flow  ( Pu  )  through  the  probe.  When  a  supersonic 
boundary  layer  is  surveyed  the  ratio  of  pitot  to  external  static  pressure 
is  usually  large  enough  (except  right  next  to  the  wall)  to  choke  the  bleed 
hole.  Hence  the  mss  flow  through  .the  probe  can  be  expressed  as: 


PuAss/>Va*  = 


m*  PpA* 

7?r 


(18) 


where  the  starred  conditions  (*)  refer  to  the  sonic  conditions  at  the 
bleed  hole  and  Hi*  is  a  Mach  number  function  tabulated  in  Ref.  39.  By 
substituting  £q.  18  into  Eq.  16  and  equating  Eqs.  16  and  17  the  following 
expression  Is  obtained: 


S*  A~  CptTr~TM*A|  s  V  (*9) 


By  combining  all  the  unknown  constants  into  one  new  constant,  B" , 


^(Tt-TM|=B"(TM-Tpl 


(20) 


Since  we  are  looking  for  small  temperature  differences,  the  variation 
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of  ./tT  can  be  neglected  leaving  the  following: 

f 

|  "p<W*BIW  (21) 

where  B  must  be  determined  from  an  experimental  calibration. 

A  probe  calibration  test  was  performed  at  Mach  6  with  the  probe 
located  in  the  free  stream  where  Tt  =  TTqw  860  R .  The  mass  flow  through 
the  probe  was  controlled  by  varying  the  tunnel  stagnation  pressure  which 
in  turn  varied  the  pitot  pressure.  The  results  of  this  calibration  are 
shown  in  Pig.  36.  These  tests  showed  that  the  measured  temperature 
decreased  (temperature  error  increased)  as  the  mass  flow  (proportional 
to  pitot  pressure)  decreased.  Furthermore ,  the  measured  variation  in  the 
temperature  error  could  be  predicted  from  the  relationship  given  in  Eq.  21. 
Also  shown  in  this  figure  is  the  temperature  error  that  would  be  predicted 
by  the  conventional  technique  of  assigning  a  constant  recovery  factor 
(  V r  =  .94)  to  the  probe.  It  can  be  seen  that  the  temperature  error 
determined  by  the  conventional  method  would  be  low  by  a  significant  amount. 
Conversely,  correcting  the  measured  probe  temperature  by  assigning  a 
constant  recovery  factor  to  the  probe  would  result  in  significant  errors 
in  the  calculated  total  temperature. 

The  method  of  correcting  the  measured  temperature  employed  Eq.  21 
with  the  constant  B  evaluated  in  the  free  stream.  The  average  probe 
temperature,  Tp  ,  was  assumed  to  be  the  laminar  recovery  temperature  on 

the  probe  head.  Thus: 
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where  Vn  i»  the  1—1  nar  recovery  factor. 
Substituting  Eq.  22  Into  Eq.  21 


Rearranging 


Pp  +  B 


U^M‘  ) 


sypP+B) 


After  nondimensionalizing  and  rearranging 


(2k) 


%  Tto 


+  B' 


’0  U-^«!  i 


Where  B;  is  the  new  constant  after  nondimensionalizing. 


(25) 


A  comparison  of  the  present  method  (Eq.  25)  with  the  conventional 
method  (  Vr  ;  const)  of  computing  the  total  temperature  variations 
through  a  typical  cooled  boundary  layer  at  Mach  6  is  shown  in  Fig.  37. 
Both  methods  correct  the  measured  temperature  at  the  edge  of  the 
boundary  layer  to  the  free-stream  total  temperature  and  both  methods 
produce  very  little  correction  near  the  wall  (correction  is  zero  when 
velocity  is  zero);  however,  between  the  wall  and  edge  of  the  boundary 
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layer  the  present  method  results  in  a  higher  total  temperature.  In 
the  integration  of  the  total  enthalpy  thickness,  *  ,  (••«  Dsta 
Reduction  Procedure)  it  is  the  temperature  difference  (l-TT/TTo )  which 
enters  into  the  calculation.  At  the  outer  edge  of  the  boundary  ***** 
the  temperature  difference  computed  by  the  present  method  is  about  half 
of  the  temperature  difference  computed  with  a  constant  recovery  factor. 
Consequently  the  value  of  *  computed  by  the  present  method  is  consider¬ 
ably  less  than  would  be  computed  by  the  assumption  of  a  constant  probe 
recovery  factor.  This  method  of  reducing  data  leads  to  a  more  conserva¬ 
tive  (lower)  estimate  of  the  ratio  of  Stanton  number  to  skin-friction 
coefficient  as  discussed  in  the  section  entitled  CORRELATION  OF  SKIN- 
FRICTION  AND  HEAT-TRANSFER  COEFFICIENTS.  The  current  evaluation  of  the 
probe  temperature  correction  is  felt  to  be  an  improvement  on  the  state- 
of-the-art  in  temperature  measurement;  however,  further  refinements  are 
still  necessary  to  improve  the  accuracy  of  temperature  measurements 
when  quantities  involving  small  temperature  differences  (such  as  ) 
are  to  be  calculated. 
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APFWDDC  II 


EEVELOFHBIT  OP  DURRBfTIAL  AID  IHTEORAL  METHODS  OP 
CORHHATEtG  SKU-PRICTIOII  AMD  HEAT-TRAHSfXR  CGEPPICUNTS 

Differential  Method 

The  following  heat-transfer  and  skin-friction  equations  relate  the 
conditions  at  the  vail  to  the  free-stream  conditions: 

S'WpfW.)8*  **-£?•  (s6) 

w 

where  the  subscript  w  refers  to  conditions  evaluated  at  the  vail.  By 
dividing  Eq.  26  by  Eq.  27  and  rearranging  the  following  form  is  obtained. 

St  _ _  ^2* 

■*T*  2(^)(T„-TW)  '  * 

Employing  the  energy  equation  in  the  following  form, 


(29) 


>>5 
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differentiating, 


du  d  u 


(30) 


and  evaluating  at  the  vail  vhere  u  =  0  ,  the  following  result 
Is  obtained: 


(31) 


By  substituting  Eg.  $1  into  Eq.  28  and  nondimensionalizing,  the  following 
form  is  obtained: 


*  (tt/tTo) 

d(u/u0) 

w 

(32) 


vhere  Pr  =  fiCp/k 

Eq.  32  can  now  be  employed  to  evaluate  St/Cf  from  the  experimental  data. 


Integral  Method 

The  following  equations  for  flat-plate  boundary  layer  growth  relate 
the  boundary  layer  momentum  thickness  and  total  enthalpy  thickness  to 
the  skin-friction  coefficient  and  Stanton  number  (heat- transfer  coeffi¬ 
cient). 

, * 

^uocpt^  (33) 

0 
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wle  *  faw**' 


(3^) 


By  dividing  Bq.  33  by  Eq.  31*  «d  rearranging,  the  following  ton  ie 
obtained. 


(35) 


If  an  isothermal  surface  (  Tw  *  const.)  is  considered,  and,  if  St/Cf  * 
const,  Eq.  35  can  be  written  in  the  following  form. 


(36) 


Then,  even  though  Cf  varies  along  the  surface,  the  integrals  cancel 
leaving  the  following  relationship: 


(37) 


For  the  case  of  a  nonisothermal  wall  which  has  a  stepwise  temperature 
distribution,  Eq.  37  can  be  written  in  the  same  manner  by  subsituting 
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an 


.all  teaperature,  r.  ,  for  T.  Ti  Ud.fl.-41v: 


Si,  ♦  fr.'V.L» 

«  ct"L,  +  cf.V 


(38) 


In  this  equation  the  .ub.crlpt  1  refer,  to  the  condition.  o.er  the 

rr:  s  t  - 

respective  portions  of  the  plate.  For  simplicity  in  evaluating  the 


rco^w*»®  —  - —  —  -  . 

data,  Cf  ,  was  assumed  equal  to  C 


*2 
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AFFUDIX  III 


1VAII2ASI0V  OF  ZHI  VALZDm  OF  US  ASSUMPTION 
FOR  THB  REDUCED  SHAPI  PARAMETER  USED  IV  US 
ANALYTICAL  nVXSTIQAZIOV 


The  method  of  analysis  used  to  predict  the  grovth  of  s  cooled 
supersonic  turbulent  boundary  layer  In  the  presence^  an 
pressure  gradient  was  based  on  the  lag-length  procedure  of  Ref.  1 
which  was  developed  for  flow  over  an  insulated  surface.  In  order  to 
modify  this  procedure  to  account  for  the  effect  of  wall  cooling  it 
was  necessary  to  evaluate  the  change  in  the  bounds  wiX 
parameter  due  to  wall  cooling.  In  order  to  provide  this  information, 
an  assumption  was  made  relating  the  shape  parameter,  6 ,  of  a  cooled 
boundary  layer  in  an  adverse  pressure  gradient  to  the  shane 

^die^tC°°i^  (adlabfttlc  vall)  tountory  layer  in  an  adver.e^r^ure 
gradient.  This  assumption  was  described  analytically  by  Eq.  8  in  the 
section  entitled  ANALYTICAL  INVESTIGATIDN.  fltf.  equ^tiL  ly  be  r^ 
written  as  follows:  ^  re 


/GCOOLED 

) 

/GCOOLED  \ 

\gadiabatic, 

/aoverse 

PRESSURE 

GRADIENT 

\gadiabatic/flat  plate 

(39) 


ine  ratio  of  the  cooled  to  adiabatic  shape  parameter  for  a  flat 
boundary  layer  was  evaluated  from  the  data  obtained  during  thTf£rt!Pii 

“  l  “d  *  P^Cr±“ZCof  the 

ratio  of  vail  temperature  to  recovery  temperature  in  Fig.  38  Both  L 

of  data  fall  along  a  single  line  with  an  scatter J i 

virtual!!,  all  the  data  lie  within  two  standard  deviations  (  2^)Tt  the 


49 


ASD  TDR  62-87 


13m  validity  of  Eq.  39  my  be  evaluated  by  comparing  the  ratio  of 
gcooleo/gadiabatic  for  *°th  the  curved-surface  data  and  the  flat-plate 
data.  If  the  data  from  both  types  of  tests  exhibit  the  saae  trends 
(when  plotted  as  a  function  of  vail  temperature  ratio),  It  aay  then  be 
concluded  that  the  assumption  of  Eq.  39  is  a  good  first  approxlaation  to 
the  physical  phenomena.  However,  none  of  the  data  were  obtained  with 
adiabatic  vail  conditions  due  to  heat  transfer  from  the  tunnel  to  the 
room.  Therefore,  before  making  the  above-mentioned  comparison,  it  was 
necessary  to  evaluate  the  reduced  shape  parameter  for  the  adiabatic  vail 
conditions . 

Tvo  different  methods  were  employed  to  evaluate  the  adiabatic  vail 
values  of  G  at  each  station  on  the  curved  surface  models.  These  tvo 
methods,  the  averaging  method  and  the  extrapolation  method,  are  compared 
in  Fig.  39  for  the  12  deg  station  on  the  Mach  3  circular-arc  model.  The 
greatest  deviation  between  these  tvo  methods  occurs  for  the  data  obtained 
at  this  station.  The  adiabatic  vail  value  obtained  by  the  averaging 
method  was  evaluated  by  employing  Eq.  39  and  inverting  the  calculation 
procedure  to  determine,  for  each  data  point,  the  value  of  Gaoiabatic 
which  is  consistent  with  the  calculation  procedure.  These  values  are 
shown  by  the  open  triangular  symbols  in  the  figures.  The  average  of 
these  values,  shown  by  the  solid  triangular  symbol,  was  used  to  determine 
the  value  of  GcoolEO/GAOiabatic  •  adiabatic  vail  value  obtained 
by  the  extrapolation  method  was  evaluated  by  fairing  a  curve  through  the 
data  and  extrapolating  this  curve  to  the  adiabatic  value  of  the  wall 
temperature  ratio.  This  value  is  shown  by  the  solid  square  symbol. 

This  method  of  obtaining  the  adiabatic  value  gives  most  weight  to  the 
data  point  at  the  highest  value  of  the  vail  temperature,  whereas  the 
averaging  method  weights  all  points  evenly. 

The  results  of  calculations  of  GCOOLEO/GAO,ABAT,C  for  the  Mach  3 
circular-arc  data  and  Mach  6  lsentroplc  surface  data  are  shown  in  Figs. 

1*0  and  hi.  In  these  figures  the  dashed  curves  cure  the  same  as  faired 
data  curves  from  Figs.  20  and  30.  As  shown  in  Fig.  40,  the  averaging 
method  results  in  several  values  for  the  reduced  shape  parameter  for 
the  adiabatic  surface  condition  at  each  station.  These  values  were 
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evenly  weighted  to  obtain  the  aean  curves  for  the  adiabatic  shape 
parameter.  These  curves  are  shown  on  the  figure  for  both  the  Mach  3 
and  Mach  6  curved  surfaces .  The  shaded  region  about  the  aean  adiabatic 
wall  curve  is  shown  to  indicate  the  maximum  accuracy  to  be  expected  on 
the  basis  of  the  2.6$  scatter  of  the  flat-plate  shape  panuneters  shown 
in  Fig.  3 6*  Moat  of  the  calculated  values  of  adiabatic  shape  parameter 
lie  within  or  very  near  the  shaded  bands.  Since  It  is  more  difficult 
to  measure  the  boundary  layer  profiles  on  the  curved  surface  models,  it 
is  to  be  expected  that  this  data  should  be  less  accurate  than  the  flat- 
plate  data.  The  results  of  the  calculations  using  the  extrapolation 
method  are  shown  in  Fig.  4l.  The  values  of  GA0,ABAT)C  obtained  in  this 
manner  always  fall  above  the  data  curves  whereas  the  values  of  GAD)ABAT,C 
obtained  by  the  averaging  method  sometime  fall  below  the  faired  curve 
through  the  uncooled  data,  see  Fig.  to. 

The  values  of  the  adiabatic  shape  parameter  at  each  station  on  the 
curved  surface  models  were  then  employed  to  obtain  the  ratio  of  GC00lE0/ 
Gadiabatic  for  each  data  point.  This  ratio  for  both  methods  is  plotted 
-in  Fig.  42  as  a  function  of  the  ratio  of  vail  temperature  to  recovery 
temperature.  The  theoretical  curve  used  in  Eq.  6  (obtained  from  the 
flat-plate  data  shown  in  Fig.  38)  is  also  shown  in  Fig.  42.  The  ratios 
computed  using  the  averaging  method  shown  in  the  upper  half  of  Fig.  42, 
scatter  about  the  theoretical  curve  with  an  r.m.s.  deviation  of  3.5$. 

The  2.8$  band  which  contained  the  flat-plate  data  of  Fig.  36  is  also 
shown  on  this  figure  for  comparison.  The  maximum  deviation  from  this 
line  is  8.6$  (or  2.5  standard  deviations)  and  occurs  for  the  Mach  3 
data  with  no  coolant  flow.  The  r.m.s.  deviation  of  3.5$  is  the  order 
of  accuracy  to  be  expected  from  the  experimental  measurements  of  boundary 
layer  profiles  in  an  adverse  pressure  gradient. 

The  ratios  computed  using  the  extrapolation  method,  shown  in  the 
lower  half  of  Fig.  42,  exhibit  a  large  deviation  from  the  theoretical 
value  at  the  low  values  of  the  temperature  ratio,  but  show  small  devia¬ 
tions  at  the  high  values  of  the  wall  temperature  ratio  due  to  the  manner 
in  which  these  points  were  weighted.  When  using  this  method  the  data 
at  a  given  station  fair  to  a  value  of  1.0  at  a  wall  temperature  ratio 
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of  1.0;  hence  the  values  of  GcocUo/^sotABATic  •* 
tore  ratios  are  forced  to  exhibit  a  relatively  lov  scatter,  For  this 
reaeon  an  r.a.e.  deviation  of  all  the  data  mold  not  represent  the 
overall  accuracy  of  Eq.  39  and  hence  <r  r.a.e.  Is  not  shorn.  In 
general,  the  ratio,  for  the  Mach  3  data  fall  hole,  the  tteoretlcal 
curve,  whereas  the  ratios  for  the  Mach  6  data  scatter  about  the 
theoretical  curve.  The  maximum  deviation  of  this  ratio  for _  the  Mach  3 
data  occur  at  the  12  and  l8  deg  stations.  Examination  of  the  data 
presented  in  Figs .  20  and  kl  shows  that  the  values  of  this  ratio  are 
s«  influenced  by  the  two  data  point,  obtained  with  no  wall 
cooling  (high  temperature  ratios).  If  the  value  of  these  two  data 

to  he  5*  lover  Uum  shown,  the  correlation 

obtained  for  the  Mach  3  data  would  be  In  ogreenent  with  the  correlation 
obtained  by  the  averaging  method. 

The  correlation  obtained  by  the  averaging  method  indicates  that 
Eq  TO  can  be  employed  to  calculate  shape  painters  within  the  limit 
of ' accurMy^of * wperimentai  data.  The  correlation  obtained  by  the 
extrapolation  method  indicates  that  Eq.  39  will  predict  the  correct 
trends  but  that  the  errors  introduced  may  be  larger  than  expected 
from  experimental  data  and  may  be  of  the  same  order  of  magnitude  as 
the  theoretical  correction  for  wall  cooling. 
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FIGURE  I  PHOTOGRAPH  OF  TUNNEL  WITH  MACH  3  NOZZLE,  COOLED 
FLAT  PLATE,  AND  SHOCK  GENERATOR  INSTALLED 
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FIGURE  2  PHOTOGRAPH  OF  TUNNEL  WITH  MACH  3  NOZZLE  AND  COOLED 
CIRCULAR  -ARC  MODEL  INSTALLED 
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FIGURE  3  PHOTOGRAPH  OF  TUNNEL  WITH  MACH  6  NOZZLE  AND  COOLED 
ISE NT  ROPIC- SURFACE  MODEL  INSTALLED 
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FIGURE  4  SKETCH  OF  HEAT-  EXCHANGER  ASSEMBLY  WITH  ISENTROPIC- 
SURFACE  MODEL  INSTALLED 
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FIGURE  7:  MACH  NUMBER  AND  WALL  TEMPERATURE  DISTRIBUTIONS  ON 
FLAT  PLATE  WITH  MACH  3  NOZZLE  INSTALLED 
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FIGURE  8:  REPRESENTATIVE  VELOCITY  AND  TOTAL  TEMPERATURE  PROFILES 
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FIGURE  II  SCHLIEREN  PHOTOGRAPH  AND  SKETCH  OF  SHOCK  -  REFLECTION 
INTERACTION  AT  MACH  3 
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FIGURE  12'  DISTRIBUTION  OF  WALL  STATIC  PRESSURE  AND  WALL  TEMPERATURE 


THROUGH  SHOCK-REFLECTION  INTERACTION  AT  MACH  3 
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FIGURE  16  SCHLIEREN  PHOTOGRAPHS  OF  FLOW  OVER  CIRCULAR -ARC 
MODEL  AT  MACH  3 
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FIGURE  I9:  TOTAL  TEMPERATURE  PROFILES  ON  CIRCULAR -ARC  MODEL 
AT  MACH  3 
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FIGURE  20-  VARIATION  OF  S*  AND  G  ALONG  CIRCULAR -ARC  MODEL 
AT  MACH  3 
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FLAT  PLATE  WITH  MACH  6  NOZZLE  INSTALLED 
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FIGURE  23:  REPRESENTATIVE  VELOCITY  AND  TOTAL  TEMPERATURE  PROFILES 
ON  FLAT  PLATE  AT  MACH  6 
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FIGURE  27  DISTRIBUTION  OF  WALL  STATIC  PRESSURE  AND  WALL  TEMPERATURE 
ON  ISENTROPIC- SURFACE  MODEL  AT  MACH  6 


FIGURE  28:  VELOCITY  PROFILES  ON  ISENTROPIC- SURFACE  MODEL  AT  MACH  6 
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FIGURE  32 ;  TYPICAL  VARIATION  OF  TOTAL  AND  STATIC  TEMPERATURES  WITH 


VELOCITY  FOR  A  COOLED  BOUNDARY  LAYER 
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FIGURE  33;  DERIVATIVE  METHOD  OF  EVALUATING  RATIO  OF  STANTON 


WALL- COOLING  PARAMETER,  (  -5-  -  * 


O  MACH  3  FLAT -PLATE  DATA 
□  MACH  6  FLAT -PLATE  DATA 


NUMBER  TO  SKIN  -  FRICTION  COEFFICIENT  FOR  FLAT- PLATE  DATA 
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AHitWiE  ?  LU  Ft 


FOR  FLAT  PLATE  HAVING  STEPWISE  TEMPERATURE  DISTRIBUTION 
THE  EFFECTIVE  WALL  TEMPERATURE,  T„'  ,  IS  DEFINED  BY  THE 
FOLLOWING  RELATION : 

w  |  *-  |  +  *-  2 

T  '  =  - ! - £ - 

L,  +  L2 


O  MACH  5  FLAT -PLATE  DATA 
□  MACH  6  FLAT-  PLATE  DATA 


WALL-COOLING  PARAMETER 


AIRFLOW 


PROBE 

STATION 


2  (TR/TT0“  Tw  /V 


=  0.805 


FIGURE  34:  INTEGRAL  METHOD  OF  EVALUATING  RATIO  OF  STANTON  NUMBER 
TO  SKIN- FRICTION  COEFFICIENT  FOR  FLAT -PLATE  DATA 
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-  SEE  FIGS. 

THAN  SHOWN  By  DATA 
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TOTAL  DEFLECTION  ANGLE 


FIGURE  35'  EFFECT  OF  ADVERSE 
STANTON  NUMBER  TO 


PRESSURE  GRADIENT  ON  RATIO  OF 
SKIN-FRICTION  COEFFICIENT 
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TfHrEHAlJ*iL  THPGR  l'&RI 
Cg*STAWT  PHQEtE  HECUVt«T 
TACTtVfl.  n  1  O-H 


RE  36:  EFFECT  of  PRESSURE  LEVEL  ON  - 

ATtlDt:  AT  M  Af.H  6 
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r«n  rnoBFrTING  PROBE  CONDUCTION  ERROR1 
ASSUMED  RELATION  FOR  CORRECTING 

Pp  (tT“tm)  5  b(tm-tp) 

WHERE  1 

B  CONSTANT  DETERMINED  FROM  MEASURED  TEMPERATURE 
WITH  MAXIMUM  FLOW  THROUGH  PROBE 

Tp  =  TEMPERATURE  of  probe  with  zero  bleed  flow 
Tt  =  Tt  a  B60  R 

I  1  ft 


5  'u 

PITOT  PRESSURE,  Pp  -IN*  H«  ABS 


* «  a  «m  inrrri  DDHRF 
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MACH  6  FLAT -PLATE  DATA 
MAXIMUM  COOLANT  FLOW  RATE 


FIGURE  37 ••  COMPARISON  OF  METHODS  FOR  EVALUATING  TOTAL 
TEMPERATURES  THROUGH  A  BOUNDARY  LAYER 


■L 
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O  MACH  3  FLAT  -PLATE  DATA 


A  MACH  6  FLAT  -PLATE  DATA 


DATA  FROM  FIGS.  10  ft  25 


ns 

o  & 

Hit  i U 

ITERS 

ON 

FLAT 

NUMBERS 

OF  3 
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O  MACH  3  CIRCULAR-ARC  DATA,  /3  =  12° 
a  GAD|ABAT|C  COMPUTED  BY  EXTRAPOLATION  METHOD 
A  GA0|ABAT|C  COMPUTED  FROM  EO  39 
▲  GAD|ABAT|C  COMPUTED  BY  AVERAGING  METHOD 

NOTE ;  THE  DATA  AT  THIS  STATION  RESULTS  IN  THE 

GREATEST  DIFFERENCE  BETWEEN  THE  TWO  METHODS 


EO  39: 


/  gCOOLED  ^ 

'  GCOOLED  j 

l  G ADI  AB ATIC  j 

CURVED 

SURFACE 

^GADIABATIcj 

PLATE 

FLAT-PLATE  VALUES  OBTAINED  FROM  FIG.  38 


FIGURE  3S-  COMPARISON  OF  EXTRAPOLATION  AND  AVERAGING  METHODS 
OF  OBTAINING  GAD|ABAT|C 
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REDUCED  SHAPE  PARAMETER 


ASD  TDR  62-87 


COOLANT 


SYMBOL 


FLOW  RATE 


JC=n 


FIGURE  40  :  ADIABATIC  SHAPE  PARAMETER  EVALUATED  BY  AVERAGING  METHOD 


DISTANCE  ALONG  SURFACE,  x  -  IN. 


SYMBOLS  DENOTE  VALUES  OF  Ga0|ABATIC  CALCULATED  FR0M  EQ'  39 
—————  MEAN  VALUES  OF 
_  _  FAIRED  DATA  FROM  FIGS.  20  8  30 


MACH  6  ISENTROPIC-  SURFACE  DATA 


MACH  3  CIRCULAR -ARC  DATA 
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VALUE  OF  GA0|ABAT|C  FROM  EXTRAPOLATION  METHOO 


FAIRED  DATA  FROM  FIGS.  20  8  30 


MACH  3  CIRCULAR -ARC  DATA 


MACH  6  ISENTROPIC- SURFACE  DATA 


DISTANCE  ALONG  SURFACE,  x  -  IN 


FIGURE  41  ADIABATIC  SHAPE  PARAMETER  EVALUATED  BY  EXTRAPOLATION 


METHOD 


93 


ASD  TDR  62-87 


■■  ■■  ■■ : 1 
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MACH  3  CIRCULAR-ARC  DATA 


MACH  6  ISENTROPIC-  SURFACE  DATA 


G AOI ABATIC  COMPUTED  BY  AVERAGING  METHOD 


rw:Q=t  ■  ICA.  n'JHUTIflK 
uHn  IS  FQ  S  BASED  On 
flat  -  pl  ate  data  or  fmj, 


rnfOtfi  ticAl  .AhiilrON 
U-5EC  'h  r-  EtAS£U  ON 

FLAT-nun:  PA r i  qf  rir, 


WALL  TEMPERATURE  RATIO,  T, 


o  D.tl— 

imimqan 

suj  snyKis;  tw;: : ; 

5  0 

CC 

O.l  0.2 

0.3  0.4 

0.5  0.6 

0.7  0 

a: 

UJ 

gadiabatic 

COMPUTED  BY 

EXTRAPOLATION 

METHOD 

FIGURE  42:  CORRELATION  OF  SHAPF  PARAMETERS  FOR  CURVED 
SURFACE  MODELS 
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